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SUMMARY 


Measurements of propeller efficiency loss due to ice formation are 
supplemented by an analysis to establish the magnitude of efficiency 
losses to be anticipated during flight in icing conditions. The measure- 
ments were made during flight in natural icing conditions; whereas the 
analysis consisted of an investigation of changes in blade-section aero- 
dynamic characteristics caused by ice formation and the resulting pro- 
peller efficiency changes . Agreement in the order of magnitude of effi- 
ciency losses to be expected is obtained between measured and analytical 
results. The results indicate that, in general, efficiency losses can 
be expected to be less than 10 percent; whereas maximum losses, which 
will be encountered only rarely, may be as high as 15 or 20 percent. 
Reported losses larger than 15 or 20 percent, based on reductions in 
airplane performance, probably are due to ice accretions on other parts 
of the airplane. 

Blade-element theory is used in the analytical treatment, and cal- 
culations are made to show the degree to which the aerodynamic charac- 
teristics of a blade section must be altered to produce various propeller 
efficiency losses. The effects of ice accretions on airfoil-6ection 
characteristics at subcritical speeds and their influence on drag- 
divergence Mach number are examined, and the attendant maximum efficiency 
losses are computed. The effect of kinetic heating on the radial extent 
of ice formation is considered, and its influence on required length of 
blade heating shoes is discussed. It is demonstrated how the efficiency 
loss resulting from an icing encounter is influenced by the decisions of 
the pilot in adjusting the engine and propeller controls. • 


INTRODUCTION 


It has long been recognized that one of the hazards of flight in 
icing conditions is the formation of ice on propellers. The presence 
of ice on propeller blades causes a decrease in the operating efficiency 
of the propeller and a corresponding decrease in the airplane performance. 
There has been much disagreement as to the actual magnitude of efficiency 
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losses experienced in icing conditions. Measurements have "been made of 
propeller efficiency losses with simulated ice formations (references 1 
and 2) and in natural icing conditions (reference 3)- The measured losses 
were low when compared with those generally expected; however, reports 
of very large apparent losses, deduced from decrease in airplane perform- 
ance, persisted. A review of the available data (references 1, 2, and 3) 
indicated that further efficiency— loss measurements supplemented by an 
analysis of the problem would be required before conclusive statements 
could be made regarding the order of magnitude of losses to 6e anticipated 
in icing conditions. The objectives of the present investigation, there- 
fore, were to provide additional measurements of propeller efficiency loss 
resulting from ice formation, and to analyze the problem of efficiency 
loss in an attempt to provide a means for applying, generally, the exist- 
ing data. Thus, additional data to determine efficiency loss were obtained 
in flight under natural icing conditions as part of a comprehensive inves- 
tigation of aircraft icing by the NACA. Another phase of the program 
(reported in reference 4) was the investigation of the meteorological fac- 
tors conducive to icing, therefore providing quantitative data on the 
severity of the icing conditions under which propeller— performance meas- 
urement's were made. 

The research was conducted by the Ames Aeronautical Laboratory. Data 
were obtained over most of the United States during the winters of 1946-47 
and 1947—48. Measurements made during the winter of 1946—47 are included 
in reference 3* The results of the 1947-48 winter season are presented in 
this report. 

Appreciation is extended to United Air Lines, Inc., the United States 
Weather Bureau, and the Air Materiel Command of the U. S. Air Force for 
aid and cooperation in the research. 


DESCRIPTION OF EQUIPMENT 


The flight tests were conducted with the twin-engine airplane shown in 
figure 1. Modifications had been made to provide thermal ice-prevention 
equipment for the wings, tail, and windshield. A description of the ther- 
mal system is given in reference 5- The right engine and propeller were 
utilized for the purposes of the research, and standard commercial elec- 
trically heated blade shoes were installed on the left propeller for ice 
protection. 

The test propeller consisted of four blades composed of double-camber 
Clark T sections. The diameter was 13.5 feet. Characteristics of the 
blade design as supplied by the manufacturer are given in figure 2. The 
symbols used in this figure and throughout the report are listed and 
defined in Appendix A. Specially built blade heating shoes, constructed 
at the Ames Laboratory, were installed for propeller heating tests, which 
were made in conjunction with the performance tests. The shoes consisted 
of several wrappings of cloth tape- in which the heater elements were 
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embedded. Two— inch— wide strips of metal were cemented over the wrappings 
along the leading edges of the blades in order to protect the shoes from 
abrasion. An enlarged spinner, built for the purpose of housing instru- 
mentation related to the heating studies, was mounted on the test propeller. 
Figure 3 shows the test propeller complete with heating shoes and spinner. 
White marking lines were painted on the blades for ease in noting ice 
locations. . 

A thrustmeter and a torquemeter were installed in the right engine 
for measuring thrust and torque absorbed by the test propeller. Both the 
thrustmeter and the torquemeter utilized pistons which were hydraulically 
restrained so that the pressure required to maintain the pistons in a 
floating condition gave a measure of thrust or torque being absorbed by 
the propeller. Pressure gages connected to the hydraulic lines provided 
means for observing and recording the values of thrust and torque. A cal- 
ibration of the thrustmeter had been made at the engine factory where the 
thrustmeter was installed. The calibration was made on a thrust dynamom- 
eter rig and gave results to an accuracy of ±1 percent through a wide range 
of propeller thrusts encompassing the normal operating range. The torque- 
meter was the standard torquemeter supplied with the engine. No special 
calibration of the torquemeter was made for the purpose of these tests. 

Since it was not feasible to calibrate the torquemeter or to recalibrate 
the thrustmeter periodically during the tests, it was decided that occa- 
sional checks of the repeatability of the thrustmeter and the torquemeter 
would serve as a satisfactory indication of their reliability. Conse- 
quently, periodic measurements were made of the efficiency of the test 
propeller in clear air. This provided a check of the repeatability of 
both' the thrustmeter and the torquemeter. The tests gave results of pro- 
peller efficiency which repeated within ±2 percent. Thus, although no 
check of the absolute values of thrust and torque was made, it is believed 
the measurements of propeller efficiency change are accurate to' within 
±U percent of propeller efficiency. 

Equipment was developed to enable photographing the test propeller 
blades while rotating so that pictures could be taken during flight in 
icing conditions without resort to feathering. The equipment consisted 
of two cameras synchronized to high-speed flash lamps and to the propeller. 
One camera was mounted to photograph the camber face and one to photo- 
graph the thrust face of the propeller blades. Figure 4 shows one of the 
cameras installed. Each camera contained two shutters, one leaf type and 
one rotary— disk type. The purpose of the shutters was to provide a means 
for eliminating as much of the background light as possible while the 
flash lamp illuminated the propeller. The disk shutter was mounted on a 
1 A— horsepower electric motor which revolved continuously during the pho- 
tographing procedure. The shutters were arranged such that the leaf 
shutter opened and allowed the aperture in the disk to pass once before 
closing again. During the time that the disk aperture was in front of 
the camera lens, the flash lamp was discharged. The duration of exposure 
was governed by the flash lamp, which gave an effective exposure time of 
about l/lO,000 second. Two flash lamps were used, one mounted beside 
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each camera, as shown in figure 5. The camera shutters, flash lamps, and 
propeller were synchronized so that pictures could he taken of any one 
selected blade with either or both cameras. Complete pictures of the pro- 
peller, consisting of both sides of all four blades, could be obtained in 
about l/2 minute. 

A resistance—" type pickup was used to indicate the propeller blade, 
angle. Calibrations showed the blade-angle indicator to be accurate with- 
in ±0.2°. 

A yaw and pitch recorder was used to record the airplane flight atti- 
tude. The yaw and pitch' recorder utilized a hemispherical differential- 
pressure— type head, as illustrated in figure 6. An electrical heater was 
installed inside the head to prevent the formation of ice which might alter 
the calibration. The accuracy of this instrument was ±l/2 . 

The meteorological instruments used to measure the liquid-^water con- 
tent, drop size, and free-air temperature of the icing clouds are described 
in detail in references 4 , 6 , and 7 * Briefly, the instruments consisted 
of the following: Rotating cylinders were used to obtain short— interval 

averages of water content and drop size. Continuous records of water con- 
tent were secured by means of a rotating-disk icing— rate meter. Free-air 
temperature was measured with a shielded thermometer consisting of a ther- 
mocouple connected to a milli voltmeter. 


TEST PROCEDURE 


Flight tests were made in clear air to obtain the datum performance 
of the test propeller, and in natural icing conditions to obtain changes 
in performance incurred through the formation of ice. The regions of oper- 
ation in natural icing conditions covered most of the northwestern, and 
part of the middlewestern. Great Plains, and eastern areas of the United 
States. Locations of the various icing encounters are listed in table I 
of reference 4 . 

During flight in icing conditions the test procedure was to collect a 
layer of ice on the propeller blades, and at the same time measure the 
liquid— water content, drop size, and free-air temperature of the icing 
clouds. Heat was applied to the wings, tail, windshield, and left pro- 
peller to keep those components free of ice. The test propeller was 
maintained in automatic pitch (constant speed) throughout the ice-build- 
up period. At the end of the period of ice build-up, the airplane generr- 
ally was climbed above the cloud layer. 

Photographs were, taken of the propeller blades while measurements 
were, made of airspeed, altitude, and attitude of the airplane; and blade 
angle, rotational speed, thrust, and torque of the test propeller. The 
purpose of making the measurements above the cloud layer was to obtain 
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stable air conditions, which were more conducive to steady readings than 
the turbulent conditions often existing in icing clouds. During each set 
of measurements the propeller blades were locked at one blade angle, and 
all flight conditions were maintained as constant as possible for a time 
sufficiently long to allow all instruments to indicate stabilized values. 
Due to the viscosity of the hydraulic oil, the pressure gages indicating 
thrust and torque took the longest time of any of the instruments to 
reach equilibrium, especially during flight at low temperatures. When 
time and other circumstances permitted, readings were taken over a range 
of advance ratio; then the blade angle was changed to a new position, 
locked, and another set of measurements obtained. In this manner, an 
effort was made to obtain data through a range of propeller operating con- 
ditions, but frequently the ice accretions would not remain on the blades 
for the lengthy period required for such measurements. Performance data 
were secured at three blade angles, 21°, 26°, and 31°> approximately, 
which covered the normal operating range of cruise and climb conditions 
for the test airplane. 


RESULTS 


A tabulation of the meteorological, flight, and propeller operating 
conditions under which propeller ice accretions were obtained is presented 
in table I. Also included is a summary of the losses in propeller effi- 
ciency obtained from each icing encounter. This col umn gives the loss in 
peak efficiency for cases where sufficient measurements were made to define 
the efficiency curve. The value was taken as the difference between the 
two curve peaks for the iced and the clean propeller. Where only point 
values of efficiency were obtained, the loss at the particular operating 
condition is given. All tabulations of conditions given in table I, with 
the exception of propeller blade angle, are those which prevailed during 
the time of accumulation of ice on the propeller. The values of blade 
angle are those maintained during the recording of performance data, after 
leaving the icing conditions. 

As will be noted in table I, small variations from the desired blade- 
angle settings of .21°, 26°, and 31° occurred due to inability to set the 
pitch exactly. In all cases, measurements of propeller efficiency were 
made at only one or two blade angles for any one particular ice formation. 

All values of liquid-water content, except for encounters 6, 8, and 10 
were measured with the rotating-disk icing-rate meter described in refer- 
ence 4. The value of water content given for encounter 6 is the average 
of three measurements taken with the rotating cylinders during the period 
of ice build-up on the propeller. The value for encounter 8 is a single 
measurement obtained with the rotating cylinders during the ice build-^up 
period; whereas that for encounter 10 is the average of four points 
obtained with the cylinders during ice build-up. 
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As indicated in table I, the test propeller was unheated except during 
two encounters. Cyclic heating was applied to the forward 10 percent of 
the blades during the ice build-up period of encounters 10 and 11. As a 
consequence of insufficient heat and inadequate chordwise coverage, 
"runback" formations collected aft of the heated areas, along with primary 
ice accumulations, but no measurable efficiency loss resulted. 

The various ice formations were grouped into different classes, 
depending on the thickness and radial extent of the formation. The clas- 
sification is based on a study of the photographs taken of the formations 
and is outlined in table II. Such a classification was established on 
the assumption that ice formations of equal thickness and extent would 
cause about the same performance change. No account was taken of the 
shape of formation since the photographs were not of such a character 
as to reveal this feature. Grouping by the size of formation alone, how- 
ever, proved to be justified, inasmuch as formations in the same class 
were found to cause approximately equal performance losses. Notation of 
the class of ice fornation is made in table I and on all the propeller 
efficiency curves. 

To establish the datum performance of the propeller with clean blades 
(no ice formations), clear-air data were secured at the three blade 
angles. The thrust— and power-coefficient curves are shown in figure 7, 
and the resulting efficiency curves are given in figure 8. Comparisons 
of the thrust— and power-coefficient curves for the propeller with the 
blades clean and with the blades iced under the conditions of table I are 
presented in figures 9(a) to 9(e). Point values of thrust and power 
coefficient for the iced propeller are listed in table III for cases 
where insufficient measurements were made to define curves. 

Curves of efficiency as a function of advance ratio for the propeller 
with ice accumulations are compared with efficiency curves for the clean 
propeller in figures 10 to 16. These curves were drawn employing values 
of thrust and power coefficients from table III and from the curves of 
figure 9 by using the formula 

1 = g 1 J (i) 

where 

tj propeller efficiency 
Crp thrust coefficient 
Cp power coefficient 


J advance ratio 
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Figure 10 presents the variation of efficiency loss for one type of 
formation at two "blade angles, and figure 11 shows the efficiency-loss 
variation for another class of formation at approximately the same blade 
angles. A comparison of efficiency losses for two classes of formation 
for the same blade angle is given in figure 12, while figure 13 compares 
losses for five formations at a higher blade angle, and figure l4 com- 
pares losses for two formations at a still higher blade angle. Figure 15 
gives the efficiency loss for a low-temperature condition in which the ice 
extended nearly to the blade tips. Figure 16 presents a plot of effi- 
ciency data for formations which accumulated during two encounters with 
propeller— blade heating. Heat was applied cyclically to the forward 
10 percent of the blade chord and, as was mentioned previously, primary 

ice formations, together with small amounts of runback, gathered back of 
the heated regions. 

It should be noted that, in figures 11, 13, and 16, singular values ' 
of efficiency are shown for the propeller in the iced condition, which 
are slightly higher than those for the clean blade. In all cases, this 
error is less than ±4-percent efficiency, which has been nbted as the 
maximum experimental error of the efficiency measurements. 


Records of yaw and pitch of the airplane showed that variations of 
2° in yaw and 1° in pitch from the clear-air conditions under which the 
datum performance curves were obtained occurred during the measurements 
of efficiency with the iced propeller. Wind-tunnel tests of the effects 
of yaw and pitch variations on propeller efficiency (references 8 and 9) 
showed that, for a propeller operating in front of a wing, changes in 
efficiency of only 1 percent resulted from yaw variations up to 5°, or 
from pitch variations up to 1°. Therefore, it is concluded that effi- 
ciency variations resulting from deviations in yaw and pitch from the 
clear -air datum conditions did not exceed 1 percent. 

Photographs of the clean propeller blades with heating shoes 
installed are shown in figure 17. These show the leading-edge metallic 
abrasion strips installed on the blades and the white— line markings 
painted on the shoes to provide a means for identifying the particular 
blade photographed and the extent of ice formations. It will be noted 
that, in some of the pictures of the thrust face of the blades, the 
leading-edge abrasion strips unfortunately give the appearance of ice 
formations and, therefore, should not be confused with the actual forma- 
tions. Pictures of the ice accretions for which efficiency data were 
obtained are presented in figures 18 to 31. All the pictures shown in 
figures 17 to 30 were taken with the photographic equipment previously 
described and with the propeller rotating. Figure 31 was obtained with 
a conventional camera and with the propeller feathered. In all cases in 
which pictures were taken with the special equipment, attempts were made 
to secure photographs of both sides of all four blades, but in many 
instances the equipment failed to function properly and only parts of the 
photographs were usable. 
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DISCUSSION 


In order to generalize the test results that have been presented, 
vith the view in mind of making possible a reasonably accurate evaluation 
of propeller efficiency losses likely to be experienced under various 
lini n g and operating conditions, the discussion that follows deals with 
(a) a consideration of those factors that influence propeller efficiency 
and their relative importance under icing conditions, (b) the changes 
that ice formations can make in these factors with consideration of the 
changes in efficiency that may be expected to result, and (c) the 
influence of operating conditions on propeller efficiency losses in icing 
conditions. 


Factors Influencing Propeller Efficiency 


The factors which affect the efficiency of a propeller may be seen 
from the following formula, which is based on blade-element theory: 


tan cp 

e tan (<P+7) 


( 2 ) 


where 


tig efficiency of blade element' 
cp angle of advance of blade element 

7 angle, the tangent of which is the 

( ten_1 lO 

An examination of equation (2) shows that, for any given operating 
condition (cp = constant), the efficiency is dependent on the blade- 
element drag-lift ratio. The primary effect of ice accretions on the aero- 
dynamic characteristics of a propeller is to increase the blade-element 
drag— lift ratio. 

Effect of radial distribution of ice accretions .- The change in 
radial distribution of thrust and torque for a propeller with a full-span 
ice formation is shown qualitatively in figure 32. The effects of a uni- 
form increase in the blade— element drag— lift ratio are illustrated in 


f tan 1 - ^ 

\ icnDx / 

blade— element drag— lift ratio 
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this figure. It is evident that the radial location of ice accumulations 
has a significant hearing on the over-all efficiency loss since the radial 
loading of thrust and torque is not uniform and the major portion of the 
loading is developed by the outer regions of the blades. An indication 
of the importance of radial location of ice accretions can be obtained from 
data of the present tests. A study of table I shows that the largest loss 
(10 percent, encounter 1) was created by a small formation which extended 
approximately 75 percent, of the blade radius; whereas a much larger, but 
shorter, formation extending about 55 percent of the distance along the 
blade caused only i)-percent loss at the same blade angle (encounter 3)o 
Thus it is evident that, when formations are limited to small radial 
extent, large accretions with correspondingly large changes in aerodynamic 
characteristics are necessary to cause large efficiency losses. An impor- 
tant factor influencing the radial extent to which ice will form on a pro- 
peller blade is the kinetic heating of the air passing over the blade. 

This effect of kinetic heating is discussed in Appendix B. 

Eff ect of changes in blade-element drag— lift ratio on propeller 
efficiency.— In order to obtain a quantitative indication of the effect 
of changes in the drag— lift ratio (represented in the symbol y in equa- 
tion (2)) on propeller efficiency, calculations were made for two hypo- 
thetical propellers for which the value of the over-all drag— lift ratio 
was adjusted to produce efficiency losses of 5, 10, 15 , and 20 percent. 

The blade-angle distribution shown in figure 2 was taken for both propel- 
lers. The two propellers consisted of conventional sections. The 
changes in drag— lift ratio for the entire blade span necessary to induce 
the assumed efficiency losses of 5, 10, 15 , and 20 percent are shown in 
figure 33 as a function of advance ratio for the two propellers operating 
at various blade angles. The envelope values of increment in drag— lift 
ratio as a function of advance ratio are shown in figure 33 for each of 
the assumed efficiency losses. It will be noted that the envelope values 
tend to decrease somewhat with increasing advance ratio. Up to advance 
ratios of about 2, the decreases in the envelope values are not great and 
the changes in drag— lift ratio required to produce the various efficiency 
losses are relatively constant. A plot of increase in drag— lift ratio for 
decrease in efficiency was made for both propellers by taking the average 
of the values of the envelope curves between advance ratios of 1 and 2. 
These data are presented in figure 3^* Propellers composed of NACA 
16 -series sections should exhibit about the same relation of change in 
drag— lift ratio to change in efficiency as shown in figure 3^» 

The data of figure 3^ apply to cases of advance ratio up to about 2. 

For advance ratios much above this value, the data of figure 3^ would no 
longer hold, due to the slope of the envelope curves shown in figure 33. 
Since the normal operating range of advance ratio for typical present- 
day transports is up to about 2, the remainder of this discussion will be 
limited to cases of advance ratio of 2 or less. 
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Effect of Ice Accretions on Airfoil-Section Characteristics 
and Attendant Efficiency Losses 


The following discussion will be divided into two parts covering 

(1) blade-element drag-lift ratio in the incompressible-flow regime, and 

(2) blade-element drag-lift ratio in the compressible-flow regime. 

Blade-element drag— lift ratio in incompressible flow.— Three investi- 
gations have been conducted which contain information bearing on this sub- 
ject. In reference 10, the effects of a mild simulated ice formation on 
the lift and drag of an NACA 0012 airfoil were studied. In this investi- 
gation, it was shown that the major change in section characteristics arose 
from an increase in drag and only a relatively small change came from 
decrease in lift. Although the accretion did not appear to be severe, the 
maximum increase in drag— lift ratio was about 100 percent, corresponding 
to a propeller efficiency loss of 8 percent. (See fig. 3^*) 

A second study, reported in reference 2, consisted of tests of a 
severe simulated ice formation on a propeller with blades composed of 
Clark T sections. The accretion extended about to the blade tips and 
was very irregular, even more so for the outer portion of the blades 
than any formation observed during the flight investigations in natural 
icing conditions conducted by the Ames Laboratory. The formation is 
believed to be at least as severe, insofar as deleterious effects on per- 
formance are concerned, as any which would occur under natural conditions. 
In reference 2, the basic data from the tests were reduced to mean lift— 
and drag-coefficient curves. These curves show the change in mean sec- 
tion lift and drag coefficient for the whole propeller and, from the 
standpoint of computing maximum propeller efficiency losses, are believed 
to represent the best data of this type available. The curves show a max- 
imum increase in drag— lift ratio of about 200 percent at a normal operat- 
ing lift coefficient (0.6), due to the presence of the simulated ice. 
Approximately l4— percent loss in peak efficiency was measured. 

The third investigation consisted of the determination of the effect 
of protuberances on the aerodynamic properties of an NACA 0012 airfoil 
(reference 11). The results of these tests, when applied to the case of 
a propeller, indicate that only very improbable' locations of ice forma- 
tions of the shape tested could cause increases in drag— lift ratio 
greater than about 150 percent. 

The results of these three tests indicate that, for propellers operat- 
ing at relatively low speeds with no effects of compressibility, the maxi- 
mum increase in blade-element drag— lift ratio possible in icing conditions 
is about 200 percent. This means, taking figure 3^ and the measured loss 
of reference 2 as a basis, the maximum efficiency loss to be expected at 
low speeds is about 15 percent. 



KACA TCT 2212 


11 


The largest measurements of efficiency loss of reference 3 are not 
ini agreement with the conclusion of a maximum loss of 15 percent. A dis- 
cussion of these measurements is given in Appendix C. 

Blade— element drag— lift ratio in compressible flow .— The only known 
high-speed data pertaining to the influence on section properties of alter- 
ations to an airfoil in the leading— edge region are given in reference 12. 
In this study, the effects of skin wrinkles on the aerodynamic character- 
istics of two NACA airfoils were investigated up to a free-stream Mach 
number of 0.73* A low-drag and a conventional section were studied. The . 
results showed that only small changes in lift resulted within the speed 
range tested. However., reductions of the order of 0,05 in drag-divergence 
Mach number were measured with both airfoils for the most adverse cases. 

In addition, the normal drag increases associated with surface irregu- 
larities prevailed over the entire speed range. From this information, 
it appears that alterations in the leading— edge contour of a propeller 
blade, of the nature that could be expected from accumulations of ice, 
could reduce the blade— element drag-divergence Mach number, thereby 
reducing the speed range for efficient operation for propellers composed 
either of low-drag (l6— series) or conventional sections. 

An illustration of the variation of drag coefficient with Mach num- 
ber such as might be expected from ice accumulations, using the skin- 
wrinkle data- as a basis, is given in figure 35; this curve is compared 
with a similar curve for a clean blade. It is apparent from these data 
that reductions in propeller performance, in addition to the performance 
losses arising from low— speed drag increases, are possible for propellers 
with a portion of the blades operating in the speed range of drag diver- 
gence. The amount of additional efficiency loss resulting from this 
cause will be a function of the length of blade affected, which is depend- 
ent on the magnitude of reduction of drag-divergence Mach number and the 
propeller advance ratio. Obviously, the greater the reduction in drag- 
divergence Mach number, the larger will be the length of blade affected. 

In order to obtain an indication of the possible magnitude of these 
increased losses, comparative calculations were made of the effects of a 
severe ice formation on the performance of a propeller operating at high 
speed ass umin g both a decrease and no change in drag-divergence Mach num- 
ber. Operating conditions were chosen to be representative of a typical 
four-engine transport airplane. The advance ratio was about 2, and it 
was established that the blade tips were operating just at the drag- 
divergence Mach number for clean blades. A decrease of 0.05 in the drag- 
divergence Mach number below that for the clean-blade case was selected 
for the iced condition. The results of the calculations showed that the 
additional decrease in efficiency due to adverse compressibility effects 
would be less than 1 percent. 

Since the tests of reference 12 on the effects of skin wrinkles 
were taken as the basis for the assumed decreases in drag-divergence 
Mach number in the calculations, the above results may not be 
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representative of the conditions prevailing for actual ice accretions . 

It is conceivable that larger reductions in drag-divergence Mach number 
may arise from the presence of ice than measured for the skin wrinkles, 
but, considering the above results, it is not likely that efficiency 
losses resulting from such causes would be greater than 5 percent. 
Therefore, it appears that for propellers operating under the most 
adverse conditions at advance ratios up to 2, with the most severe full- 
span ice accretions possible, efficiency losses would not exceed 20 per- 
cent. 


Normally expected efficiency losses.- In the previous sections, it 
was concluded that the maximum efficiency loss to be expected at low 
speeds is about 15 percent and, for propellers operating in the critical 
speed range, the maximum loss would not exceed 20 percent. For the major- 
ity of cases of propeller icing, the losses would be less than these 
values since they are representative of extreme conditions in which the 
entire length of blade is adversely affected; whereas, in a large number 
of instances, ice forms on only part of the blade and usually in a much 
less detrimental configuration.. This is in general agreement with the 
measured values given in table I, which displays relatively low losses, 
with a maximum value of 10 percent. Also, a review of the measurements 
made during the flights of reference 3 showed that, for about 90 percent 
of the icing encounters, efficiency losses were less than 10 percent. 

Thus it appears that, in the vast majority of instances of propeller 
icing, efficiency losses can be expected to be less than 10 percent. 

Mention was made previously that large efficiency losses have been 
reported, based on decreases in airplane performance. In some cases these 
reported losses have been as high as 30 percent. Such large values are 
in conflict with both calculated and measured losses, and it is believed 
that accumulations of ice on other components of the airplanes caused 
serious increases in drag with corresponding decreases in airspeed, 
which were falsely attributed to propeller icing. Measurements reported 
in reference 3 showed that increases in airplane parasite drag of as 
much as 80 percent are possible when all components, with the exception 
of the propellers, have accumulated ice. 


. . Influence of Operating Conditions on Propeller 
Efficiency Losses in Icing Conditions 

When a propeller accumulates ice, the resulting changes in propeller 
performance are reflected in corresponding changes in airplane perform- 
ance. Depending on how the pilot reacts to the changes, various operating 
conditions can be established which may have an increased detrimental 
effect on the propeller performance. This aspect of propeller icing can 
best be explained by reference to a plot of efficiency as a function of 
advance ratio (fig. 36). In this figure, curve A represents the envelope- 
efficiency curve of a clean propeller; whereas curve B represents the 
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equivalent curve for a case where ice accretions would cause an efficiency- 
loss of a fixed amount. The fact that the two curves would he approxi- 
mately parallel is substantiated by the data of reference 1. If an air- 
plane operating at point C on curve A encounters icing and the pilot does 
not adjust his controls, the airspeed will decrease and the airplane will 
operate at some point such as D, assuming the propeller has a constant 
revolution— per-minute pitch control. Upon noting this loss of airspeed, 
the pilot may take two courses of action. He may try to return to a 
point near the peak of the envelope, such as point E, by adjusting 
.engine speed and power. If the airframe has accumulated ice, resulting 
in higher airplane drag, this would probably call for a reduction in pro- 
peller speed and an increase in engine power. The pilot may try to shed 
some of the ice on the propeller by increasing propeller speed. The com- 
bined action of centrifugal force and kinetic heating (see Appendix B) 
resulting from an increase in propeller rotational speed is often effec- 
tive in reducing the extent of the ice accumulation. If he is unsuccess- 
ful in this respect, the propeller will continue to operate on curve B at 
some point such as F. He has thus lost additional efficiency with respect 
to the possibility of operating at E, shown as in figure 36, due to 

the curvature of the envelope curves. On the other hand, if some of the 
ice is removed, he would operate at some point G, which approaches 
curve A as a limit. Then, by resuming the initial engine operating con- 
ditions, it may be possible to operate at or near the original point C. 

It should be noted that the foregoing discussion applies to an unpro- 
tected or an inadequately protected propeller. If the entire airplane is 
properly protected, there will be no need for corrective action by the 
pilot since the propeller will continue to operate near point C. 

Thus it appears that in operation of unprotected or inadequately 
protected propellers in icing conditions, periodic attempts should be 
made to throw, off the accretions by increasing propeller speed. If the 
ice cannot be removed, it is desirable to operate at or near the peak of 
the efficiency-envelope curve for the iced propeller. This peak is 
probably located at about the same value of advance ratio as the peak 
.for clean blades. 


CONCLUSIONS 


As the result of measurements of propeller performance loss due to 
ice formation, made during flight in natural icing conditions and sup- 
plemented by an analysis of the factors contributing to performance loss 
(in which all values apply for propeller advance ratios up to about 2), 

the following conclusions are reached: 

/ * 

1. During the vast majority of icing encounters with unprotected 
propellers on present-day transports, propeller efficiency losses can be 
expected to be less than 10 percent. 
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2. The maximum loss to he anticipated for propellers free of adverse 
compressibility effects is about 15 percent; whereas, for propellers oper- 
ating above the critical speed, the maximum loss to be expected is about 
20 percent. 

3. .Reported propeller efficiency losses larger than about 15 or 
20 percent, based on reductions in airplane performance, probably are 
actually due to icing of other components. 

4. The magnitude of efficiency loss resulting from an icing encounter 
can be influenced by the decisions of the pilot in adjusting the engine 
and propeller controls. 


Ames Aeronautical laboratory. 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., June 15, 1950. 
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APPENDIX A 
SYMBOLS 


chord of blade element, feet 

specific heat of air at constant pressure, Btu per pound 
blade— element drag coefficient 



propeller diameter, feet 
drag of blade element, pounds 

acceleration of gravity, equal to 32.2 feet per second, second 
blade section maximum thickness, feet 

mechanical equivalent of heat, equal to 778 foot-pounds per Btu 
propeller advance ratio 

lift of blade element, pounds 
free— stream Mach number 



propeller rotational speed, revolutions per second 
input power to propeller, foot-pounds per second 
torque of propeller, foot-pounds 

radius to blade element, feet 

/ ' 

propeller tip radius, feet 
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t temperature, °F 
T thrust of propeller, pounds 

U speed of free stream with respect to propeller-blade section, feet 
per second 

V airplane true airspeed, feet per second 
x fraction of tip radius 



P blade angle, degrees 

P Q blade angle at 0.75 radius station, degrees 
7 angle between lift force and resultant force ( tan 


H propeller efficiency 


Ct'N 


T| e efficiency of blade element 


tan .tp 


Do \ 

T> 


degrees 


L tan (9+7 ) J 

p mass density of air, slugs per ’cubic foot 

cp angle of advance of propeller— blade element at radius x 


/ 




tan 


— l 




JtriDx/ 


, degrees 
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APPENDIX B 

CONSIDERATIONS OF THE RELATION OF KINETIC 
HEATING TO PROPELLER ICING 


Kinetic— Temperature Rise Experienced. By a Propeller 
During Flight in Icing Conditions 


In clear air the temperature rise at the stagnation point of a blade 
section is given by 


At = 


U 2 

2gjc p 


(3) 


This equation, which is derived from a balance of the thermal and mechani- 
cal energies involved, defines the air— temperature rise, and, provided 
there is no thermal conduction in the propeller blade, the surface temper- 
ature also will experience the same increase. If, as is true in the prac- 
tical case, there is conduction in the blade, the surface temperature rise 
will be less than that given by equation (3) by an amount dependent on the 
equivalent conductivity of the blade, its shape, and other factors. 

In a cloud composed. of water drops, the kinetic— temperature rise of 
the blade surface is reduced due to cooling by evaporation of the drops 
impinging on the surface. Equations for computing the kinetic— temperature 
increase with evaporation prevalent are presented in reference 13. One 
difficulty involved in the solution of the equations of reference 13 is 
that they must be solved by trial, a laborious procedure. The stagnation- 
temperature rise in a cloud may be obtained with considerably less effort 
by using a pseudo-adiabatic diagram, such as shown in figure 37- The 
resulting value will be the same as that obtained using the method of ref- 
erence 13. The procedure for computing this temperature rise is outlined 
in f igure 37 ' as follows : 

Establish point A at the pressure altitude and ambient-air temperature 
for which the calculation is to be made. Calculate At^^ using 
equation (3), add this value to the air temperature, and follow down a 
dry adiabatic line to the temperature line obtained from the resulting 
sum (point B). This establishes the pressure in terms of altitude at 
the stagnation point. Point C is then determined by the intersection 
of the stagnation— pressure line and the line drawn through point A 
parallel to the pseudo-adiabatic lines, and this point is the wet 
stagnation temperature. It will be noted that the wet kinetic— temper- 
ature rise is considerably less than that obtained in dry air. 
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There are two limitations to the above method for computing stagnation- 
temperature rise in wet air. First, no account is taken of the heat removed 
from the stagnation region by water not evaporated which flows aft during 
flight in clouds of high water content. This factor likewise is not con- 
sidered in the analysis of reference 13 . The second limitation is that heat 
conduction in the blade is neglected. There is no simple way of including 
this factor in the calculations. A fairly complete analysis and discussion 
of this aspect of kinetic heating is included in reference 13 . In refer- 
ence 13 , comparisons are made between the effects of kinetic heating on a 
perfectly conducting blade and a completely nonconducting blade in a cloud. 
Only conduction in a chordwise direction is considered. It is shown that 
the effect of conduction is to reduce the stagnation temperature from that 
for a nonconducting blade. 

In order to obtain a comparison of kinetic— temperature rise calculated 
by the above methods with temperature rises prevailing during actual flight 
in icing conditions, values obtained during the present tests were plotted 
in figure 38 , together with curves computed using equation (3) and a pseudo- 
adiabatic diagram. The atmospheric conditions assumed for the curve com- 
puted using the pseudo-adiabatic diagram are 12,000 feet pressure altitude 
and 5° F free-air temperature, which are the average of the conditions for 
the flight values of kinetic— temperature rise. The values obtained during 
the tests were computed using propeller photographs from icing encounters 
in which the radial extent of the ice formations apparently was limited by 
kinetic heating and in which there was no obvious breakoff. By assuming 
that the end of the formation was at freezing temperature and noting the 
radial extent from the photographs, it was possible to plot the kinetic- 
temperature rise above ambient as a function of the section velocity at 
the end of the formation. For further comparison, experimental data from 
reference 14, which were acquired in a manner similar to that employed for 
the above points, also were plotted in this figure. 

In the case of the values obtained during the present investigation, 
it was noted that the points displayed considerable scatter. A study of. 
the data revealed that the points fell into three groups, depending on the 
atmospheric conditions prevailing at the time the photographs were taken 
which established the radial extent of the ice formations. The grouping 
is given in the following table and the data points of figure 38 were cor- 
respondingly divided: 


Group 

Atmospheric condition existing 
at time of photograph 

1 

Continuous icing 

O 

Icing clouds of very low liquid- 
water content or broken clouds with 
patches of clear air between 

3 

Clear air or snow 
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Group 1 covers the points for the condition in which sufficient liquid 
water would he present to insure continuous evaporation at the blade lead- 
ing edge. The points encompassed by group 3 should exhibit close to the 
full kinetic— temperature rise for clear air in the absence of evaporative 
cooling from water on the surface. In group 2, periodic exposure of the 
blade to water clouds, then to clear air, or exposure to very small amounts 
of water, would cause a smaller degree of evaporation than created under the 
conditions of group 1. Thus, the points of group 2 would be expected to 
fall between the points of groups 1 and 3- Such actually proved to be the 
case. In figure 38 , the points falling in group 1 established a fairly 
well-defined curve, those in group 3 were distributed around the curve 
determined by equation (3), and the points in group 2 lay in between the 
two curves. With one exception, the data points from reference l4 followed 
the same general pattern, covering the same range as the other points. The 
nonconforming point obtained at a blade-section velocity of 500 miles per 
hour is believed to be in error. 

It will be noted that the curve for continuous icing (group 1) is 
below that established by the pseudo-adiabatic diagram. This could be due 
to a combination of two causes. First, some of the kinetic heat developed 
may have been removed by water flowing back from the leading-edge region. 
Second, thermal conduction in the blade would tend to lower the experimental 
values. Neither of these factors is considered in the procedure for deter- 
mining the calculated curve. If conduction were the predominant cause of 
the lower experimental values, this effect also would be prevalent for the 
clear-air points (group 3)> but such is not the case. The material used 
in the construction of the test heating shoes appears to have provided 
fairly good insulation. It is concluded, therefore, that removal of heat 
from the leading-edge region. by water flowing aft can have considerable 
effect in reducing the kinetic—' temperature rise of a propeller blade in 
icing conditions. 


Beneficial Effect of Kinetic Heating in Preventing 
Premature Compressibility Effects 


With the outer portion of a propeller operating at high velocity, it 
is possible that the deleterious action of ice accretions in reducing the 
blade-section drag-divergence Mach number would be prevented through the 
beneficial effects of kinetic heating. Thus, if ice is prevented from 
forming on the blades in the critical regions, no adverse results from 
this cause can occur. The extent to which the beneficial effects will be 
manifested is a function of the propeller speed and the free-air tempera- 
ture. A plot of the air temperature down to which protection from kinetic 
heating will be obtained as a function of the blade— section Mach number is 
given in figure 39* Values of kinetic— temperature rise from the continuous- 
icing curve (group 1) of figure 38 were used in the construction of 
figure 39- 



20 


MCA TFT 2212 


The data of figure 39 niay be used to determine the ambient tempera- 
ture above which adverse compressibility effects due to the presence of 
ice will be prevented by kinetic heating. If a drag-divergence Mach 
number of 0.75 is assumed (typical of the -outer sections of present-day 
propellers), kinetic heating will prevent compressibility effects result- 
ing from ice accretions down to an ambient temperature of about 2° F. 


Length of Blade Heating Shoes as Influenced 
by Kinetic Heating 


In most cases at the present time, propeller -ice protection is 
obtained through the use of rubber heating shoes cemented to the external 
surface along the blade leading edges. The shoes are subject to erosion 
by the abrasive action of foreign matter in the air. This is. especially 
true for the outer regions of the blade where velocities are high. The 
problem is so critical from a maintenance standpoint that most airline 
operators have resorted to using relatively short heating shoes to keep 
abrasion at a practical minimum. A typical extent of shoe currently in 
use on a 13— foot-diameter propeller is to the 60-percent— radius station. 
For such a propeller operating at 1000 rpm and 300 miles— per— hour true 
airspeed, complete protection could be anticipated down to free-air 
temperatures of about 15° F, due to the combined action of the heating 
shoes and kinetic heating. Below this temperature ice could form on the 
blades in the region limited on the inner end by the blade shoes and on 
the outer end by kinetic heating. At very low temperatures, of course, 
the ice would extend from the outer end of the blade shoes clear to the 
tips. Under these circumstances, it can be shown, using blade-element 
theory, that, with an ice formation over the unprotected portion of the 
blade of sufficient severity to cause an over-all efficiency loss of 15 
percent for a full-span accretion, a loss in efficiency of about 10 
percent could result. Thus it is seen that partial— span blade shoes can 
reduce considerably, but by no means entirely, propeller efficiency 
losses resulting from the accretioh of ice. In a good many instances, 
air temperatures will be sufficiently high to prevent any ice formation 
beyond the blade shoe, and only the negligibly small efficiency losses 
resulting from periodic formations on the blade shoes through cyclic 
heating will be incurred. 
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APPENDIX C 

CONSIDERATION OF PREVIOUS MEASUREMENTS OF LARGE PROPELLER 
EFFICIENCY LOSSES IN ICING CONDITIONS 


Three measurements of propeller efficiency loss presented in refer- 
ence 3 produced values of 15 percent or greater. The air temperatures 
prevailing at the time of these measurements were relatively highj thuB 
the action of kinetic heating (see Appendix B) would have limi ted the 
radial extent of the ice accretions. Under such circumstances, the pre- 
sence of adverse compressibility effects is extremely unlikely, and the 
maximum efficiency loss to he expected would he about 15 percent, as 
concluded in the body of the present report. With only partial-span 
formations, it would be expected that the measured losses should have 
been somewhat less than 15 percent. In view of this discrepancy between 
the measured losses of reference 3 and the maximum loss as concluded in 
thiB report, it is of interest to investigate the changes in prope ll er 
section properties required to produce the measured losses. Accordingly, 
calculations were made to establish the magnitude of aerodynamic changes 
necessary to cause the three losses. 

The three highest losses recorded in reference 3 are 15, 17, and 19 
percent. The air temperature existing at the time of the measurements 
was 25° F, 21° F, and 18° F for the losses of 15, 17, and 19 percent, 
respectively. From this information, the radial extent of the formations 
prevailing for these temperatures was computed, using the curve of 
kinetic-temperature rise for continuous icing (group 1) of figure 38. 

The airplane operating conditions given in reference 3 were used in the 
computations. Results of this calculation established the outer limit 
of the ice accretion at 33 percent of the blade radius in the case of the 
15-percent Iobb, at 60-percent radius for the 17-percent loss, and at 

69-percent radiuB for the 19-percent loss. 

\ 

The changes in section drag— lift ratio which would be required for 
the iced portions 'of the blades to cause the measured efficiency losses 
were then computed, using clean-blade thrust— and torque-distribution • 
curves which were considered representative. For the 15- and 19-percent 
efficiency— loss cases, these curves were calculated by the method of 
reference 15 > and, for the 1 7— percent loss case, they were approximated 
from data of reference -l6. Results of the calculations showed that over 
2000— percent increase in drag-lift ratio would have been required for 
the iced part of the blades to cause an efficiency loss of 15 percent 
under the existing conditions. Similarly, increases in drag-lift ratio 
of 900 and 800 percent for the iced regions were calculated to have been 
necessary to create the losses of 17 percent and .19 percent, respectively. 

In view of the fact that an increase in drag-lift ratio of about 
200 percent was concluded previously in this report to be the maximum 
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possible in icing conditions, it is difficult to reconcile measurements 
producing 4 to 10 times this change. Possibly the formations encountered 
during these particular tests of reference 3 were as deleterious as the 
measurements indicate, but the validity of these three values is seriously 
questioned. It is not believed that losses of such a magnitude can be 
experienced at the relatively high air temperatures of the tests. During 
the experiments of the present report it was learned that erroneous meas- 
urements of efficiency are very likely if considerable care is not exer- 
cised in allowing conditions to stabilize before taking readings. This 
is sometimes very difficult to accomplish, in view of the turbulence often 
associated with icing conditions, and it is believed that the measurements 
of the three efficiency losses under discussion are in error. Validity of 
the results of reference 2, from which the maximum value of 200-percent 
increase in drag— lift ratio was derived, is not disputed since the tests 
were conducted in a wind tunnel in which stable conditions generally exist 
also, there is no reason to believe the data were subject to systematic 
error, and a large number of observations were made, which would tend to 
reduce random errors. 
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TABLE II 

ICE FORMATION CLASSIFICATION SYSTEM USED THROUGHOUT REPORT 


Class 

Outer limit of ice 
formation 

(percent tip radius) 

Class 

Thickness of ice at 
leading edge in chord- 
wise direction (in.) 

I 

26 to kl 

A 

Less than l/4 

II 

42 to 54 

B 

l/U to l/2 

III 

55 to 68 

c 

l/2 to 1 

IV 

69 to 8l 

D 

Over 1 

• 

V 

82 to 95 

— 

— 
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TABLE III 


EXPERIMENTAL VALUES OF THRUST AND POWER COEFFICIENT FOR SEVERAL 
ICING ENCOUNTERS IN WHICH ONLY SINGULAR VALUES 
OF PROPELLER EFFICIENCY WERE MEASURED 


Encounter 

number 

Advance 

ratio, 

v/nD 

Thrust 

coeffi- 

cient, 

C T 

Power 

coeffi- 

cient, 

C P 

Blade 

angle, 

(deg) 

Data 

figure 

number 

Photo 

figure 

number 

4 

1.07 

0.066 

0.078 

25 

.8 

11 

2 

2 

j 

3 


1.07 

.064 

.077 

25 

.8 





i 


1.16 

.104 

.159 

31.8 





1 


1.17 

.101 

.152 

31.8 





i 

f 






V 

\ 

/ 



• 95 

.074 

.084 

25 

.8 

13 

25 



.98 

.074 

.084 









1.03 

.071 

.086 




. 


• 



1.03 

.070 

.083 


' 





V 




\ 

/ 

V 

V 

8 

1.07 

.073 

.088 

25.8 

13 

26 



1.07 

. .073 

.088 









1.07 

.073 

.089 



. 






1.08 

.073 

.088 









1.08 

.072 

.088 









1.09 

.073 

.089 







V 




\ 


V 

V 

9 ' 

.90 

.104 

.115 

21.5 

12 

27 


.92 

.109 

.121 

21.5 

V 

V 

10 

1.13 

.096 

.126 

26.5 

1.6 

28, 

29 



1.14 

.090 

.118 







\ 

( 

1.15 

.096 

.127 

■| 

\ 

f 


\ 

/ 

n 

1.28 

.102 

.149 

29.5 

16 

80 

vi 

1.29 

.102 

.149 

29.5 

V 

\ 

/ 

12 

• 98 

.115 

.139 

28.0 

15 

31 



1.00 

.111 

.l4o 









1.04 

.102 

.128 

j 






V 



i 

\ 

i 

\ 

/ 

: Y 1 
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Figure 1.— The test airplane used for the flights in natural icing conditions to investigate propeller- 
efficiency loss. 
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Figure 2. - Characteristics of the te 
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Figure 3 .— Test propeller .shoving the specially constructed blade heating 
shoes and the spinner. 
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Figure 4 .— Camera installation for photographing the thrust face of the 
test propeller blades. 



Figure 5.- Propeller camera-flash-lamp installations. 
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Figure 6.- Hemispherical differential-pressure type yaw and pitch head. 
Airspeed total— pressure heads are mounted on either side of the yaw 
and pitch head. 
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Figure 8.- Efficiency curves for the test propeller with blade heating 
shoes and with no ice formations . 
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Figure 9. - Continued. 






Figure 9.- Continued. 
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Figure 10. - Comparison of efficiencies for dean and iced propeller 
at two blade angles for ice formation type HD. Thrust- and 
power - coefficient data from figure 9(a). 
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Advance ratio , J 


Figure i i.~ Comparison of efficiencies for dean and iced 
propeller at two blade angles for ice formation type 
ID + i/eRA. Encounter 4. Thrust - and power- coefficient 
data from table JR. 
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Figure 12.- Comparison of efficiencies for dean 
and iced propeller for two ice formations. 

Thrust- and power- coefficient data at Q s 20.7° 
from figure 9(b), at Q = 2i.5° from table HI . 
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Clean 

Iced, form, type JTC * 1 / 4 YB, enc. 6 

Iced, form, type HD, enc. 345 

— Iced, form, type HB+JHA, enc. t 


O Iced, form, type I A, enc. 7 
□ Iced, form, type HA, enc. 8 



.8 .9 t.O l.i 1.2 t.3 

Advance ratio , J 


Figure 13.- Comparison of efficiencies for dean and 
iced propeller for five ice formations at one 
blade angle. /3 o - 25.8° . Thrust- and power- coef- 
ficient data from figures 9 (a) , 9 (c), and 9 (d), and 
table HI. 
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Figure 14.- Comparison of efficiencies for dean 
and iced propeller for two ice formations at 
one blade angle. )9= 3 i.8°. Thrust - and power - 
coefficient data from figures 9(a) and 9(e). 
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Advance ratio, J 


Figure 15 - Comparison of efficiencies for dean and 
iced propeller for ice formation type TEB. fi o = 28.0°. 
Encounter 12. Thrust - dnd power -coefficient data 
from table JR. 
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Advance ratio, J 


Figure 16 - Comparison of propeller efficiencies 
for blades dean and with smalt ice formations 
resulting from inadequate cyclic heating . Thrust - 
and power -coefficient data from table HI . 
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(b) Thrust face of blade. 

Figure 17.— Photographs of clean propeller blades. 
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Figure 18.— Ice formation of encounter 1. Camber face. Peak efficiency 
loss, 10 percent. 
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Figure 19-— Ice formation of encounter 2. Camber face. Peak efficiency- 
loss, 1 to 2 percent. 
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Figure 20.- Ice formation of encounter 3 . Camber face 
loss, k percent at p Q = 25. 8°; 6 percent at 0 O = 


Peak efficiency 

31 . 8 °. 






Figure 21.— Ice formation of encounter 3a. Camber face. Peak efficiency 
loss, 0 to 1 percent. 
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(a) Camber face. 


Figure 22.— Ice formation of encounter 4. 
0 O = 2 5.8°; 8 percent at 0 O = 31 . 8 °. 


Efficiency loss, 0 percent at 
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(b) Thrust face. 


Figure 22.— Concluded. 
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(a) Camber face. 

Figure 23.— Ice formation of encounter Peak efficiency loss, 
4 percent. 
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(b) Thrust face. 


Figure 23.— Concluded. 
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(a) Camber face. 

i 

Figure 24.— ice formation of encounter 6. Peak efficiency loss, 
4 percent. 






(t>) Thrust face. 
Figure 24.— Concluded. 
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(a) Camber face. 



(b) Thrust face. 

Figure 25.— Ice formation of encounter 7* Efficiency lo6s, 1 to 
3 percent. 
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(b) Thrust face. 

Figure 26.- Ice formation of encounter 8. Efficiency loss, 0 percent. 
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(a) Camber face. 

Figure 27.- Ice formation of encounter 9. Efficiency loss, 1 to 
4 percent. 
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(h) Thrust face. 
Figure 27.— Concluded. 
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Figure 28.— Ice formation of encounter 10. Thrust face. Efficiency loss 
0 percent. 
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Figure 29.— Ice formation of encounter 10. Thrust face. Efficiency loss 
0 percent. 



Figure 30.— Ice formation of encounter 11. Thrust face. Efficiency loss 
0 percent. 
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Figure 31 *— Ice formation of encounter 12. Camber face. Efficiency loss, 
5 percent. 
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(b) Torque . 


Figure 32- Qualitative indication of changes 
in spanwise distribution of thrust and 
torque for a propeller with a full- span 
ice formation. 
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Advance ratio, J 
(b)NACA 0012. 

Figure 33.- Calculated curves of change in section drag -lift ratio for entire blade 
span required to produce various efficiency losses for two propellers composed 
of different airfoil sections. 




Increase in drag -lift ratio, 


MCA TTf 2212 


91 



Change in propeller efficiency, A ip 


Figure 34.- The calculated variation of the maximum 
decrease in propeller efficiency with increase in 
over -all blade - element drag -lift ratio. Data not 
applicable for advance ratios greater than 2. 
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ft ‘AOUBtOfJJd J9U9dOJd 


Figure 36,- Comparison of various propeller operating possibilities in conditions 
of icing. 



Pressure altitude , ft 
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Dry adiabat 

Pseudo- adiabat 



Figure 37.- Pseudo- adiabatic diagram showing procedure for 
computing kinetic-temperature rise at stagnation point 
in a cloud (wet air). 
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Nonconducting blades, calculated from equation (3) 

Nonconducting blades, calculated using pseudo - 

adiabatic diagram (fig. 37) 12,000 feet altitude , 

5° F air temperature 

— O — Continuous icing (group!) 

□ Low liquid -water content or alternate 
dear air and cloud (group 2) 

O Clear air or snow (group 3) 


} Data this 

report 


A Flight data, ref. 14 




Figure 38. - Comparison of calculated and measured values of 
kinetic -temperature rise at stagnation line of a propeller 
blade as a function of blade -section velocity. 
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Jo * 9JfUDJ9dUJ94 JjD - 99JJ 




NACA-Langley - 10-26-50 - 1150 


Figure 39- Variation with blade- section Mach number of free - 
air temperature down to which ice formation will be 
prevented by kinetic heating. 
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Pll-X 

012 

013 

• * 01 n 



021 

$ 22 "^ 

023 

* * 02 n 



031 

032 

833 -X . 

• • 03n 

= 0 

( 11 ) 

• • 0 
Pm 

#. • • 

P n 2 

* * * * 

Pn3 ‘ 

• • • • 




or 


where 



( Ha) 


= B ° A_ 

ihq 5 4 to 2 


The dynamical matrix is 


M • 


• Coupled hendlng-torslon vibration , - The torsional and tending 
deflections due to coupled bending-torsion vibrations of a canti- 
lever beam are given by (see appendix B): 


where 
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absolute magnitude of projection of distance from elastic axis 
to center of gravity or. perpendicular to bending direction 
for root section 


rgQ radius of gyration about elastic axis at root section 

Equations ( 2 ) and ( 8 ) define and Pjj- The quanti 

ties and ere given by 


ru = T1 K jk - Of- 1 ) s k M' Jk . + ^ s r m 

k= 1 k 


r=k+l 


j 


where 

p k 

P 1k s . V " (k* 1 ) 2 + J(k-1) 2 f.(z)dz 
J* ^k-i L 2 2 J J 

ok 

Q jk s J |(k-l) 2 z + |(k-l ) 3 fj(z)dz 

- 


Sjf is the ratio of average static mass unbalance of 
k^k interval to static mass unbalance at the root section. 

The Station Numbers P . ^ and are listed in tables I 

to VIII with the other Station Numbers. The determinantal equa- 
tion becomes 



NACA TN 2185 


11 


ra ir 

-X Pa 12 

. . ra ln 

cTrn 

€rn 2 . 

• € ^^In 

ra 21 

Fa 22 - X 

• • Fa 2n 

€r ?21 

*^722 • 

* € ^2n 

• • 
ra nl 

• • • • 
ra n2 

• • 

nn 

• • • 

€r ^nl 

,r 'n2 ' 

• eI \n 

5 11 

8 12 

• * 8 In 

Pll"* 

Pl2 . 

• • Pin 

8 21 

8 22 

* * 8 2n 

^21 

^22"^ * ' 

’ • ^2n 

• • 4 

8 nl 

8 n2 

* * 8 nn 

• • • t 

p nl 

» • • • • 

Pn2 * • 

• * Pnn"^ 


or 


1 -hu]| “ 0 ( 14a ) 

vhere jr^jJ is the dynamical matrix and I is the identity matrix. 

The first n roots of equation (14) will give the first 
n coupled frequencies. 


APPLICATIONS AND RESULTS 

In applying the previously discussed method, it is necessary 
to determine for a given team the elements ajj, p.^, and 

of the dynamical matrices. These quantities will depend on 
the physical properties of , the beam and on the number of stations 
chosen. If the physical properties of the beam are known, the quan- 
tities a^j, PiJ> 7ij> an( i 8 ij can 8 ® directly calculated from 
equations (2), (8), and (13). The numbers Mj k , Nj k , P Jk , Q Jk , 
M 'jk> N ' jk> 8114 Q’jk a PP®aring in these equations depend 

on the number of stations n that are used and can be read directly 
from tables I to VIII for any given number of stations up to eight. 
Once these quantities have been calculated, equations (6), (11), or 
(14) can be solved for the frequencies by any method desired. The 
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matrix-iteration method used herein is simple and rapid and 
requires no particular computing skill. As will he indicated, 
however, the accuracy of equations (6), (11), and (14) is such 
that relatively few stations need he used, in which case it may 
he convenient to expand the determinants and to solve the result- 
ant low-order algebraic equation. 

In order to illustrate the accuracy, this method was applied 
to torsional vibrations, bending vibrations, and coupled vibra- 
tions of a uniform cantilever beam. The exact theoretical values 
for torsional vibrations and bending vibrations of uniform canti- 
levers are well known. The exact theoretical values for the coupled 
bending-torsion vibration of a uniform beam were calculated (ap- 
pendix D) . A comparison was then made between the values obtained 
by the method presented and the exact theoretical values. The 
number of stations used was 1, 2, and 3 (n = 1, n = 2, and n *= 3). 
The comparisons are summarized in table IX. 

Torsional vibration . - For the case of a uniform beam, 

= 1^ s 1 and equation (2) becomes 


i n 

a ij =2Z! N jk - (fc- 1 ) + y~ Mj r (15) 

k=l l_ r=k+l 

The values of and are given in tables I to Till. 

The table to be used depends on the choice of the number of 
stations. 

Let n = 1 

'’•an = N n 

From table I, = 5/l2 

* 5/l2 


± I* JO M 2 9l 

1 12 C 0 


or 
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0 ) 


2 12 C 0 


5 T ?2 
I 0 t 


2.400 




0) = 1.549 



The exact theoretical value for the first torsional frequency is 

(0= 1.571 

The error is -1.4 percent when only one station is used. 



The mode shape obtained by the method of Station Functions 
agrees .well with the theoretical mode shape, as is shown in fig- 
ure 2(a). 

Let n m 2, then by equation (15) and table II 


a ll “ W 11 + M 12 


_§_ + _5_ = 57 
15 12 = 60 


a 12 " N 21 + M 22 = 


31 + 29 _ 57 

240 48 120 


a 21 ■ N 11 + K 12 


_8_ 8 16 
15 15 “ 15 


a, 


22 


N, 


21 


N. 


22 


31 | 239 _ 13 
240 240 = 15 


The determinantal equation then becomes 


57 _57 

60“ A 120 


16 13 

15 15 “ 


0 
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which gives 

= 1.6214 
Xg = 0.1953 

Therefore 

(0 X = 1.571 
(0 2 = 4.526 

The exact theoretical values are 

CO x = 1.571 



The percentage error of the first two modes, for only two sta- 
tions, is found to he 0 and. -4 percent. 

The mode shapes are shown in figures 2(h) and 2(c). Agree- 
ment of the first mode with the exact theoretical shape is excel- 
lent; the second mode agrees fairly well. 

Let n = 3, then hy equation (15) and table III 

a il “ N 11 + M 12 + M 13 “ °- 945833 

a-12 = Ngi + Mgg ■*" Mgj * 0.958333 

a !3 - N 31 + “32 + “33 " °‘ 520834 

Ogi = N n + H lg + 2M 13 




St 


1.033333 
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Ogg « ^21 ^ + 2 M 23 = 1.883333 

a 23 - N 31 + N 32 + 2M 33 " ^O 11113 
a 31 = N 11 + N 12 + N 13 * 1*012500 
a 32 “ N 21 + N 22 + N 23 * 2 * 025 0 00 


a 33 = W 31 + N 32 + N 33 “ 1* 3875 °1 


The determinantal equation is ' 

0. 945833 -X 0.958333 

1.033333 1.883333-X 

1.012500 2.025000 

The solutions are 


0.520834 
1.011113 
1.387 501- X 


0 


Xj = 3.6474 
X g ■* 0 . 4093 
X 3 = 0.1599 

Therefore 



I 
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The exact theoretical values are 



The percentage errors of the first three modes, calculated by use 
of three stations, are found to be 0, -0.5, and -4.5 percent, 
respectively. 

The mode shapes are shown in figures 2(d) to 2(f). The first 
two modes agree very veil with the theoretical shapes; agreement of 
the third mode is fair. 

This procedure can be carried out as shown for any number of 
stations desired. 

Bending vibrations . - For a uniform beam, Bjj. = m^ = 1 and 
equation (8) becomes 



and from table I 

Pr.-Ji-.31 59 

J " L 630 1008 720 
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Therefore, from equation (7) 

CO » 3.493 

The exact theoretical value is 

CO *» 3.516' 


JO 

PO * 4 

w 

lv 4 


The percentage error for just one station is found to he 
-0.65 percent. 

The mode shape is shown in figure 3(a) and is seen to agree 
very well with the theoretically exact shape. 

Let n = 2; then by equation (16) and. table II 


Pii = p, n-^ii + F , i2 


-¥ 


*12 = 0.422745 


»12 * P ’21 - «'21 + ? M ’22 -^'22 = °' 29592S 
e 21 - 2P ’ll * 2 P 'l2 - «'n - «’l2 4 H 'l 2 


1.145167 


022 « 2P? 21 + 2p, 22 " V21 ‘ ^22 + | N ’22 -| M '22 = 0.905530 


The characteristic equation is 

|o.422745-A 


The roots are 


1.145167 


0.295925 
0.905530- A 


* 1.2943 
A g “ 0.0339 




The exact theoretical values 
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The characteristic equation is 

0. 270604 -X 1.009943 0.487441 

0.648170 3.266250-X 1.689891 = 0 

0.985135 5.822852 3.204301-X 

The roots are 

X l = 6.5521 
X 2 = 0.1667 
X 3 =0.0223 


Therefore 


The exact values are 




The percentage error for three stations is found to be 0, 0, and 
-2.4 percent, respectively. The modes are plotted in figures 3(d) 
to 3(f). The first two modes are seen to agree very well with the 
theoretical mode shape; agreement of the third mode is fair. 
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Coupled bending-torsion vibrations . - A uniform team vith 
the following constants was chosen: 

7 = = 38.56 

€ «= 0.8 


r 

e r 


193.2 


n L 


241.5 


The values of and 3^ 
the same as given before for 
because = B k «= C^ a m k ■ 


are obtained as previously and are 
n = 1, n b 2, and n = 3. Also, 
Ijj. = 1, equation (13) becomes 




ij 


k=l 


N' 


Jk 


(MJ "** *£&“•* 


'u-h*** 


VL. 

Q Jk + Z _ 


r=k+l 


<l-ta|) B Jr ) 3 . 


M 


Let n = 1, then the determinant is 


Jr 


ra ll‘ * 


tTy 


5 11 

The roots are 


11 


Pit* 


0.002 156 -X 

0.111111 


X ! = 0.0837 
X 2 — 0.0005 


0) 


1 “ 3.46^ 


B, 


»2 = 44.7' 


mol- 

B 0 

/V 


0.001196 
0. 081944 -X 
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The procedure for calculating the exact theoretical values is 
derived in appendix D. The exact values are 


(0 ^ 3 .49 


B 


0 




0 ) 2 







The percentage error for the first mode, calculated "by use of one 
station, is -0.9 percent. 


Let n «= 2, then the determinant is 


ran-* 

Fa. _ 
12 

err u 

€ r ^12 

r a 

21 

Ta - X 

22 

tTy 

'21 

er 7 

'22 

8 11 

5 12 

Pll“X 

Pl2 

6 21 

8 22 

^21 

022’^ 


Substituting the known values and solving for X gives for the 
first two roots 


Xl = 1.3197 
X 2 « 0.0412 
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and the frequencies "become 



The percentage errors for two stations are: 

-0.3 percent for the first mode 
-4.4 percent for the second mode 

This procedure can "be carried out for any number of stations 
desired. For three stations, the frequencies obtained are 



, The percentage errors are 

-0.3 percent for the first mode 
0 percent for the second mode 
-1.8 percent for the third mode 

The results obtained by the method presented are seen to agree 
very well with the exact theoretical values. 
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These results are summarized in table IX, where a comparison 
is also made with the results obtained for uncoupled bending and 
torsional vibrations by use of influence coefficients with weighted 
matrices (reference 12). The values using weighted matrices were 
taken from table I of reference 12. It can be seen that for a 
given number of stations, the results* obtained by the method pre- 
sented herein are considerably better than those obtained by using 
influence coefficients with weighted matrices. In general, it is 
indicated that for a uniform cantilever beam using n stations 
along the beams, the first n - 1 frequencies and modes are in 
excellent agreement with exact theoretical values and even the n th 
mode is given within the accuracy with which the physical properties 
of the material are known. For a tapered beam, more stations may 
be required, depending on the amount of taper. The number of stations 
required to give satisfactory accuracy is listed in table X. A com- 
parison is also made by using weighted influence coefficients; the 
values are taken from table II of reference 12. 


The first vibrational frequency is also given approximately 
by equation (C2) appendix C) when coupling exists between bending 
and torsion; it is plotted in figure 4. In order to check these 
curves, the exact solution was obtained (appendix D) for the ratio 

(C0 t /0 b ) 2 equal to 4 and was plotted on the same figure. The val- 
ues given by equation (C2) are seen to be in excellent agreement 
with the theoretically exact values. 

The effect of the coupling between bending and torsion is to 
reduce the first natural frequency below that which would exist if 
there were no coupling. This effect is shown in figure 4, wherein 
the value of 0 is always less than 1. This decrease in the first 
natural frequency due to coupling is, however, relatively unimportant 

in the practical range of '((O^/o^ ) ^ > 4 and e < q.75. 


SUMMARY OF RESULTS 

A method based on the use of Station Functions is presented 
for calculating uncoupled and coupled bending -tors ion modes and 
frequencies of arbitrary continuous cantilever beams . The results 
of calculations made by this method indicated that by the use of 
Station Functions derived herein, at least n - 1 modes and fre- 
quencies can be obtained accurately by using just n stations 
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along the beam if the beam is uniform. For a tapered beam, more 
stations may be required, depending on the amount of taper. The 
amount of computational labor is markedly less than for other 
methods. The use of Station Numbers tabulated herein further re- 
duces the amount of calculation necessary. It is shown that the 
effect of coupling between bending and torsion is to reduce the 
first natural frequency to a value below that which it would have 
if there were no coupling. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio. October 18, 1949. 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 

&ij coefficient in equation for Station Function 

in torsion 

B bending stiffness of beam, function of z 

B 0 bending stiffness at root section of beam 

ratio of average bending stiffness of 

k^ 11 interval to bending stiffness of 
root section 

coefficient in equation for Station Function 
in bending 

C torsional stiffness of beam, function of z 

Cq torsional stiffness of root section of beam 

Cjr ratio of average torsional stiffness of 

*fcll 

k interval to torsional stiffness at 
root section 


c -l> c 2> c 3 - constants defined in appendix B 

fj( z ) Station Function in torsion for sta- 

tion (defined in text) 

gj( z ) Station Function in bending for sta- 

, tion (defined in text) 

I mass moment of inertia per unit length of 

beam about elastic axis, function of z, 
except where otherwise defined 

*0 mass moment of inertia per unit length of 

beam about elastic axis at root section 
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Ik 


i, j,k,n 

j, k,r 


ratio of average mass moment of inertia 
per unit length of k til interval to mass 
moment of inertia per unit length at root 
section 

station indices 
summation indices 


l length of beam 

M jk> N jk> p jk> Qjk> Station Numbers (defined in text); function 

^'jk> W jk> I”jk> Q'jk °I indices j, k, and n 

m mass per unit length of beam, function of z 

mass per unit length of beam at root section 

% ratio of average mass per unit length of 

k th interval to mass per unit length at 
root section 


n 


1b( z ) 


r go 


number of stations along beam 

bending loading function on beam 

torsional loading function on beam 

absolute magnitude of projection of distance 
from elastic axis to center of gravity on 
perpendicular to bending direction 

radius of gyration about elastic axis at 
root section 


r 0 absolute magnitude Qf projection of dis- 

tance from elastic axis to center of 
gravity on perpendicular to bending 
direction for root section 

S static mass unbalance, function of z, mr 

static mass unbalance at root section, mQTQ 
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S 


k 


x 


7 

*i 

z 

a i y p ij’ 7 iy 
b ±y Vj 


ratio of average static mass unbalance at 
k section to static mass unbalance at . 
root section 

distance from root of beam, except where 
otherwise defined 

bending deflection, function of z 
bending deflection at i th station 
dimensionless distance along beam, x/6 

elements of dynamical matrix defined in 
text 


r 

7 

5 

c 

e 

e i 

* 

Q 

03 

(0 b 


“t 


_i_ £o 
8 2 c 0 mo 

uncoupled frequency ratio (O3 b /oo b ) 2 

length of interval along beam between two 
stations 

p 

coupling coefficient ( r o/ r g o) 

torsional deflection, function of 1 z 

torsional deflection at i th station 

root of frequency equation or character- 
istic root of dynamical matrix 

2 

frequency ratio (to /oo^ ) 

frequency of vibration . . 

frequency of uncoupled fundamental bending 
mode 

frequency of uncoupled fundamental torsional 
mode 

second derivative of deflection with respect 
to time 
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APPENDIX B 


STATION FUNCTIONS AND DETERMINANT AL EQUATIONS 
Torsional Vibrations 

A schematic diagram of a cantilever beam divided into n inter- 
vals of length 5 is shown in figure 1. The Station Functions for 
the torsional vibrations of such a beam must satisfy the follow- 
ing conditions: 

at 


z = 0 

f i (o) = 0 

(Bl) 

z b n 

f ' j(n) b o 

(B2) 

z . i 

f i(i) » 1. 

(B3) 

z = i 

- 0 j/i 

(B4) 


where f'(z) denotes the derivative with respect to z. 

Equations (Bl) and (B2) represent the boundary conditions that 
must be satisfied by a cantilever beam vibrating in torsion; equa- 
tions (B3 ) and (B4) represent the further conditions imposed upon’ 
the Station Functions. These conditions will be satisfied by a 
function of the type 


f l (z) - a li z * a 2i z2 + • • • + a (n t l)i z<n+1> 


(B5) 


where the coefficients a^ must satisfy the following simul- 
taneous equations obtained from conditions (B2), (B3),and (B4): 

0 « a^ + 2na 21 + 3n 2 a 31 + . . . + (n+l)n n a( n+ i) 1 (B2a) 

1 = ia n > i 2 a 21 + i3a 31 + . . . + i(» +D a (a+1)l (B3a) 

0 - J a xi + J 2 a 21 + J 3 a 3i + . . . + j( n+1 ) ± j/i 


(B4a) 
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The coefficients a^ can "be obtained by solving equations (B2a) 
to (B4a) and the functions f^(z) determined for each station. 
Equation (B5), however, can also be written in the following form; 


f i(z) = 


jji ^ Z "*^ z ( z “ c l) 


n (i-J) i(i" c i) 


(B5a) 


where ^ represents the product for all values of j except 

J * i. The function in equation (B5a) obviously satisfies condi- 
tions (Bl), (B3), and (B4), because it has zeros at all points 
specified by equation (B4), it equals 1 at the point specified by 
equation (B3), and it equals zero at the point specified by equa- 
tion (Bl). In order to satisfy condition (B2), the constant c^ 
is determined by substitution of equation (B5a) into equation (B2), 

c-^ = n for l/n 


c^ = n [1 + 




for i = n 


Equation (B5) can be obtained from equation (B5a) by carrying 
out the indicated multiplications. The complete deflection function 
is then given by 

e(z) » f 1 {z)e 1 + f 2 (z)e 2 + . . . + f n (z)e n 

n (B6) 

J=1 

The continuous loading function q t (z) can now be written as 
q t (z) a Iu?e{z) a Id) 2 i f V(*)0J (B7) 


A continuous loading function, which is a function of the 
deflections at the reference stations, has thus been obtained. 



30 


NACA TN 2185 


Bending Vibrations 

The Station Functions for the bending vibrations of the beam 
shown in figure 1 must satisfy the following conditions: 


z « 0 g^O) = 0 (B8) 

z = 0 g'i(O) = 0 ( B9 ) 

z = n g’ ^(n) = 0 (BIO) 

z = n g' ' ^(n) = 0 (Bll) 

z = i gi (i) =1 ' (B12) 

z = J gi(j) = 0 j/i (B13) 


where g'(z), g”(z), and g"*(z) denote the first, second, and 

third derivative, respectively, of g(z) with respect to z. 

Equations (B8) to (Bll) represent the boundary conditions that 
must be satisfied by a cantilever beam vibrating in bending and 
equations (B12) and (B13) represent the additional conditions 
imposed upon the Station Functions. 

These conditions will be satisfied by functions of the type 

. gi (z) = b 2i z 2 + b 31 z 3 + ...+• ^ (n+3)l z (n+3) (B14) 

where the coefficients b^j must satisfy the following equations 
obtained from conditions (BIO) to (B13) : . 

0 = 2bg^ + 6n bj^ + . . . . . . + (n+3) (n+2)n^ n+ -'-) ^( n+ 3 )i 

(BlOa) 

0 . 6D S1 4 24nl 41 ♦ . . . + (n+3) (n t 2) (n.l)n“b (n+3) j (Blla) 

1 * b 21 + l 3 b 31 + ■ • • + i ( "* 3) t( n+3)l (B12a) 

0 * J 2 * J 3 » 3 i + • • • + J (n+3) b (n+3) i J^ 1 


(B13a) 
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The coefficients can therefore he obtained from equations (BlOa) * 
to (B13a) and the functions gi(z) determined for each station i. 
Equation (B14) can, however, be written in the following form; 

w, (z-j) z 2 (z 2 +c_z+c,) 

8i(0 - = -,- ■■■ (B14a) 

(i-J) i (i +c 2 i+c 3 ) 

, n 

where ^ represents the product for all values of j except 
«J = i. The function in equation (B14a) obviously satisfies con- 
ditions (B8), (B9), (B12), and (B13), because it has zeros at all 
points specified by conditions (B8), (B9), and (B13) and equals 1 
at the point specified by equation (B12). In order to satisfy ' 
conditions (BIO) and (Bll), the constants c 2 and c 3 are deter- 
mined by substitution of equation (B14a) into equations (BIO) 
and (Bll). The general forms for c 2 and c 3 are, however, com- 
plicated and it is easier to obtain the numerical values of these 
constants for each specific case. Equation (B14) can then be 
obtained from equation (B14a) by carrying out the indicated multi- 
plications. The complete deflection function is then given by 

n 

x( z ) = XZ 8<( z ) 7 . (B15) 

J=1 J J 

The continuous bending loading function q b (z) can now be written 
as 

q b (z) = mofy(z) = mof ^ g.(z)y (B16) 

j=l J J 


Coupled Bending-Torsion Vibrations 


The Station Functions for the coupled bending torsion vibra- 
tions are the same as previously given for the bending vibrations 
and the torsion vibrations. The loading functions, however, are 
given as follows (see reference 7): 

, O^U) = ld> 2 0(.z) + SO) 2 y(z) 


(0 


£ 


Tfj(z)0j + Sgj(z) 


(B17) 
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and 


2 2 

q- (z) = So) e(z) + moo y(z) 


n " 

£ sf 

j=l L 


2_, | s fj( z ) 0 j + “s^ 2 ) yj 


(B18) 


Determinantal Equations and Dynamical Matrices 


Once the Station Functions and the corresponding loading 
functions have "been determined, the deflections at the reference 
stations can he obtained in terms of the loading function. A 
homogeneous equation in the reference-station deflections for each 
station is thereby obtained. The determinant of the coefficients 
of the resultant set of homogeneous equations can be set equal 
to zero; the determinantal frequency equation is thus derived. 

The deflections at the reference stations are obtained by the well-' 
known equations for obtaining influence coefficients. 


Torsion. - The deflection at the station i due to the con- 


tinuous loading q^(z) on the beam is given by 


0 ^ = 8 




q t ( z ) 


n n 


dz 


i. dz '+ 5 2 


^0 C 


<M Z ) 


dz 

i dz (B19) 

^0 C 


If C is assumed to have a constant value for each interval, 
then these integrals may be written as the sum of integrals over 
each section. Equation (B19) then becomes 


1 Co 4 ti Ck 


IJk-1 


z qfc(z) dz + 


nk 


u k-l 


(1-k) q t (z) dz + 


O n 




q t (z) dz 


(B20) 


By substituting the relation 

q t (z) = to 2 
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and by assuming a constant value for I for each interval and 
changing the summation order 



Ik 


k 

p r 

r=k+l '-'r-l 


zf (z)dz - (k-1) I k f,(z)dz + 

1 J vk.x J 


f j(z)dz 


(B21) 


Let 


n k 




Uk-1 

^k 

Jk-1 


z V z) dz= V 


f j ( z ) dz= Mj k 


(B22) 


J 


Then 


£ »* £ i 4 v 


0 .3=1 


(B23) 


where 




T k M jk + XZ! I r M jr 
r=k+l 


(B24) 


If C* - I k - 1 (constant cross section), then 


“U-C 


ktrl L 


n 

N Jk - (k-1) M 1k + YZ M 


Jk 


r=k+l 




Let 


X = 


'0 


_ 2 2 


(B25) 


(B26) 
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Then 


X& 1 = a i j 0 J 


(B23a) 


and the characteristic equation is 


N • 


M = 0 


(B27) 


where I is the identity matrix. 


Bending . - The deflection at the station i due to the con- 
tinuous loading q b (z) on the beam will be given by 

pi oZ 

4 (z-z )(i-z ) 

y = 5 q. (z) dz- dz + 

1 ' uq “ Oq B 1 


4 , , (z-z^U-z-^ 

5 <l b (z) dz dz 

i A D \Jn B i 


(B28) 


If B is assumed to have a constant value for each interval, 
these integrals may be written as the sum of integrals over each 
interval. Equation (B28) then becomes 


f ph 

r i = — XU — I 1 4- - (k-l)z + 

B 0 fci B k [ ^k-i L 2 

Vk . 1 [r~ I' 11 - 1 ) 2 * + J* 1 '- 1 ) 3 ] 


)z + i(k-l) q b (z)dz - 


q fc (z) dz + 


' n p n r t 3~| 

|z - |(2k-l)J q b (z) dz + |(2k-l)z - j q b (z) dz 

(B29) 


By substituting the relation 
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lb 


(z) = 0) 2 


m 


£ 

. 1=1 


gj (*)>j 


(B30) 


and "by assuming a constant average value for m in each interval 
and changing the summation order 


y< = 


0) 2 m () 5 4: n 


E». 


1 B 0 ft iJ J 


(B31) 


where 


H) =4^ "r 


(l- k 4> H ' Jr + 


^ ■) "V. 


(B32) 


vk 


p» = 

J k o k .i 


Q, J k ”vJk-l 


\ - (k-1) z + |(k-l) 2 




gj (z) &z 


r xk 


N 1 = 
w Jk- 


Wk-1 J 
nk 


Y - |(k-l) 2 e + i(k-l) 3 


Z g, (z) dz 


gj(z) dz 


M'., = 


g. (z) dz 


(B33) 


J 


^k-1 • J 

For a uniform beam, = B^ « 1 and equation (B32) becomes 
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P 


-£ t 


n 

u i ip, jk-«’ 3 k + m 
ld k=l V 0 0 r=k+l 


|(i-k+g-) N' 


Jr + 


Is 


^-(k-l) 3 (2k-l) 


M' 


Jr. 


(B32a) 


Let 


B, 


X a 


C 0 2 6 4 m„ 


(B34) 


then the characteristic equation becomes 


W- 


XI 


0 


where I is the identity matrix and jp ^ , j 
In expanded form, equation (B35) becomes d 


(B35) 

is the dynamical matrix. 


Pn-x 

P 12 . * 

* * P ln 



p 2i 

• 1 • • 

^22“^ * 

• • • • • 

* * P 2n 

• i • 

= 0 

(B35a) 

Pnl 

P n2 

* « 




where X is a latent root of the matrix 


M- 


Coupled bending-torsion vibrations . - The deflections at the 
station i are given as before by equations (B19) and (B28). The 
loading functions q t and q^ are changed, as follows: 


q t (z) = to 2 |jE 0(z) + S y(z)J 
q t (z) = (o 2 [s 0(z) + m y(z)"| 


(B36) 
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If these two equations are substituted into equations (B19) 
and (B28) and the integrations are performed as previously, the 
following relation is obtained: 


B, 




T- S K y i . 




y< W 2 m 0 5 4 ^ / y 

= “57“ ( &i J 0 J + — 


(B37) 


0/ J 


L 0 0 J=1 

where and are given in equations (B24) and (B32) and 


s ^0_ 

r e r 


r ^ 

& 2 Cof^o 


ij 


,Eff 


n 


s k N' Jk -■ X^S r M' Jr 


» * 

^ i) « j 


(i-k^) N Jr + 


(B38) 


f j (z) dz 




>k 

r , 1 j - If 11 - 1 ) 22 + |< k -i) 3 ] fj urn* 


where 
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the determinantal equation therefore is 


XI - 



vhere [^ij] is ‘ b ^ ie dynamical matrix, 
as indicated in -equation (B37). Ft). ,1 
matrix. L JJ 


the elements of vhich are 
is seen to be a 2n x 2n 
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APPENDIX C 


QUADRATIC FORMULA FOR FIRST COUPLED MODE 

If only the first vibrational mode is desired, it is possible 
to obtain this mode approximately by coupling together the funda- 
mental uncoupled bending mode vith the fundamental uncoupled tor- 
sional mode to obtain a simple quadratic equation for the first 
coupled frequency. This equation is valid when the coupling coef- 
ficient c is constant along the beam. The differential equa- 
tions obtained by coupling the fundamental uncoupled torsional mode 
with the fundamental uncoupled bending mode 

2 

my + S0 + ma^ y * 0 
Sy + 10 + I O) t 2 0 = 0 

where 

m mass per unit length of beam, function 

S static mass unbalance, function of z 

I mass moment of inertia about elastic axis, function of z 

CO-jj frequency of uncoupled fundamental bending mode 

0 >t frequency of uncoupled fundamental torsional mode 

• • denotes differentiation twice with respect to time 

These equations lead to a quadratic equation in the frequency 
ratio 0 , whose solution for the lowest frequency, provided € is 
constant along the beam, is 

0 . iu_ 

<V 2(1-0 

where 

ft , frequency ratio, (o)/o> b ) 2 
y uncoupled frequency ratio, 


*7(1-6) 


1 - 


(i-r)2 


( 02 ) 


are: 


of z 


(Cl) 
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c coupling coefficient, (r/r_) 2 

O 

This quadratic has heen plotted in figure 4 for values of e 
ranging from 0 to 1 and values of y = (w^/to^) 2 from 1 to 100. 
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APPENDIX D 


EXACT SOLUTION FOP COUPLED BENDING-TOPSION VIBRATIONS OF 
UNIFORM CANTILEVER BEAM 

The differential equations for the equilibrium of an element 
of a beam vibrating in coupled bending -tors ion vibrations can be 
put in the following dimensionless form: 


where 


d%2 ml 


ml ■ 




dx 


* 2 Y 2 

dx 2 


r Y 1 + <0° Y 2 

B B 


IV 


■ - € 


Yi - 


Il‘ 


■CO Xr 


J 


Yq= y/r 
y 2 = 0 

= distance from root 


(Dl) 


now 


€ = (r/r g ) 2 


0)- 


2 = fii 

5 ml^ 


V = =5 4; 

II" 


where 


c 4 s 12 *36 


c 5 = '2.467 

Equations (Dl) become 
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where 


Let 


Then 


ft 4 *! 

dx 4 

d 2 Y 2 
dx 2 = 


~\ 

c 4 fl(Y 1+ Y 2 ) 

) 


- e 




Q 5 (0)/(0 b ) 2 

7 = (^t/ w h) 2 


«1 

dx 


*3 





dYc 

if = c 4 ^(Yl+Y 2 ) 



'5Q 


< £ VV 




(D2) 


(D3) 


Equations (D3) can he written as the single matrix equation 
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0 0 1 0 0 0 

0 0 0 0 0 1 

.0 0 0 1 0 0 

0 0 0 0 1 0 

c 4 Q c 4 Q 0 0 0 0 

-€05^ -c 5 Q 0 0 0 0 


where Y and A are the matrices indicated. 

The solution to the matrix equation (D4) is given by 

Y = e^ Y Q 

where Yq is a column of arbitrary constants. 

From the boundary conditions 


at x = 0 


Y 3 - 0 


X *r 1 


Y= - Y c = 0 


Y n = Y(0) * 


0 

0 

0 

y 4 (0) 

y 5 (0) 

y 6 (o) 


If then Q is an element of the matrizant e A , the boundary 

1 J 

conditions give 
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044 045 0 46 

054 055 056 =0 (D6), 

064 065 066 

Equation (D6) is the frequency equation. It has an infinite num- 
ber of roots for 0). 

In order to determine the elements 0 4 j, must be evalu- 

ated. Use vill be made of Sylvester's theorem (reference 13). 

0 

The X matrix of A is 

-X 0 10 0 0 

0 -X 0 0 0 1 

0 0 -X 1 6 0 

0 0 0 -X I 0 

c 4 0 c 4 0 0 0 -X 0 

c^3 c 5 0 0 0 0 -X 

7 € 7 

The characteristic equation A(X) * 0 is 

X 6 + ^ X 4 - c 4 0 X 2 - ( 1-c ) c 4 c 5 = 0 (D7) 

2 

Equation (D7 ) is a cubic equation in X . Let the roots be 

' ”^1> ^ 2 * ^3* “^3 

Then by the confluent form of Sylvester's theorem 
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From equations (D9) and (DIO), the elements Q,, are seen to 
be given by J 

3 4 °^ 


*44 


Z — i — p _ o cosh 

1=1 ifx - X J 5 


^ 45 = " 


Q 


4 + \ 2 
1 _y A 1 


sinh 


i=l X.* XI n 2 4 2i 

Vi <*1 ~h > 

V c.QX , 2 

n 46 = - 2 _j T ? o- sinh X. 

46 fel A i n (X ± 2 -X J 2 ) 1 




C/C 


c 4 nx 1 2 + — l - (i-c) 


n 2 


« 54 =-E 


-x o p sinh X. 

i=i A i_n (X 4 -xp 1 




Q 55 = fi 44 




^c R b - / , ■ " - — “ — - — posh X. 

56 ui n (\ 2 \ 1 

j/i <*i -*j > 


Q 64 


3 -€c —X. 

— 57 i 


- / — — sinh X. 

1=1 'Si <** v> 


0 


65 


.-r 


“ €c 5 ? A i 


^66 “ 


I^ A iII (X^-X 2 ) 
J^i D 

X i 4 " C 4 Q 


cosh Xj 


i=i n a 2 -x 2 


Si <v -V) 


(Dll) 


cosh 
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The value of the determinant in equation (D6) must be - plotted 
against the frequency; the value of the frequency for which this 
determinant becomes zero is thereby obtained. This procedure 
involves first solving the cubic equation (D7) for each assumed 
value of frequency parameter and then calculating the elements of the 
determinant from equations (Dll). The process is evidently long, 
and laborious. 
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TABLE I STATION NUMBERS 
n = 1 




1 

M 

1 

2 

3 

N 


5 

12 

P 


3 

20 

Q 


7 

180 

M' 


2 

5 

N' 


13 

45 

P» 


71 


TO 

Q’ 


31 

1008 


TABLE II - STATION NUMBERS 
n = 2 




i 

2 

M 

l 

11 

12 

5 

15 

N 


8 

8 


15 

15 

p 


0.183333 

0.025000 

Q 


.046032 

.029365 

M* 


.536364 

.627273 

N* 


.367100 

.851948 

P* 


.137933 

.057955 

Q* 


.036616 

.069733 

M 

2 

.13 

29 



48 

48 

N 


- 31 

239 


240 

2^0 

P 


-0.037500 

0.143750 

Q 


-.008135 

.181448 

M' 


-.060795 

.448674 

N» 


-.034875 

.758685 

P' 


-.011252 

.118462 

Q' 


-.002614 

.150415 
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TABLE III - STATION NUMBERS 
n = 3 



X 

1 

2 

3 

M 

1 

0.950000 

0.450000' 

-0.050000 

N 


.545833 

.587500 

-.120833 

P 


.186310 

.032143 

-.005357 

Q 


.046577 

.038244 

-.011756 

M» 


.596268 

.533205 

-.097426 

N' 


.399646 

.708205 

-.239994 

P» 


.148013 

.042560 

-.012798 

Q» 


.038884 

.050843 

-.028318 

M 

2 

-0.525000 

0.725000 

0.475000 

N 


-.241667 

1.175000 

1.091667 

P 


-.068452 

.160714 

.031548 

Q 


-.014583 

.202083 

.068750 

M» 


-.149356 

.602896 

.625418 

N' 


-.033406 

.994860 

1.475153 

P' 


-.026378 

.143948 

.057937 

Q’ 


-.006034 

.181698 

.127659 

M 

3 

0.235185 

-0.153704 

0.568519 

N 


.106019 

-.231944 

1.513426 

P 


.029563 

-.023677 

.139749 

Q 


.006222 

-.028963 

.316408 

M» 


.040630 

-.072928 

.445812 

N» 


.022325 

-.111744 

1.200133 

P» 


.006972 

-.012081 

.118007 

. Q’ 


.001579 

-.014830 

.267865 
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TABLE IV - STATION NUMBERS 
n = 4 




1 

2 

3 

4 

M 

1 

1.022222 

0.429630 

-0.051852 

0.022222 

N 


.576455 

.557937 

-.127249 

.076455 

P 


.194478 

.029597 

-.006581 

.002612 

Q 


.048240 

.035167 

-.014547 

.008359 

M* 


.623188 

. 511882 

-.082891 

.042276 

N» 


.413738 

.676680 

-.203719 

.146954 

P» 


.152256 

.039616 

-.010651 

.005795 

Q' 


.039818 

.047267 

-.023551 

.018630 

M 

2 

-0.647917 

0.747917 

0.518750 

-0.085417 

N 


-.292857 

1.207143 

1.207143 

-.292857 

P 


-.081920 

.163021 

.041295 

-.009598 

Q 


-.017295 

.204828 

.090642 

-.030688 

M* 


-.211987 

.667412 

.544025 

-.112648 

N» 


-.116662 

1.091462 

1.269193 

-.390585 

P» 


-.036502 

.153469 

.044508 

-.015000 

Q* 


-.008281 

.193310 

.097745 

-.048203 

M 

3 

0.522222 

-0.255556 

0.633333 

0.522222 

N 


.229365 

-.381746 

1.673810 

1.729365 

P 


.062798 

-.037401 

.148512 

.037202 

Q 


.013040 

-.045624 

.335791 

.118397 

M» 


.122052 

-.166738 

.582158 

.643846 

N» 


.065879 

-.252823 

1.545802 

2.164827 

P» 


.020304 

-.026235 

.140822 

.060554 

Q* 


.004551 

-.032114 

.318707 

.194016 

M 

4 

-0.221701 

0.094850 

-0. 105961 

0.543924 

N 


-.096544 

.140724 

-.267841 

1.997206 

P 


-.026267 

.013391 

-.017322 

.136803 

Q 


-.005428 

.016301 

-.038574 

.446729 

M» 


-.035456 

.042628 

-.064723 

.438962 

17* 


-.019023 

.064205 

-.164169 

1.622066 

P' 


-.005836 

.006481 

-.010869 

.117037 

Q’ j 

-.001303 ' 

.007917 

-.024222 

.382752 
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1*097991 

.608222 

.202887 

.049943 

.649902 

.427616 

.156411 

.040729 


-0.839550 

-.373049 

-.103119 

-.021583 

-.255330 

-.139170 

-.043239 

-.009758 


0.762798 

.329315 

.089079 

.018334 

.197103 

.105126 

.032115 

.007151 


-0.548214 

-.233780 

-.062665 

-.012809 

-.117990 

-.062435 

-.018958 

-.004201 


0.214238 

.090928 

.024289 

.004952 

.033722 

.017789 

.005389 

.001192 


2 

3 

4 

5 

0.408755 

-0.040898 

0.019866 

-0.013120 

.527493 

-.100112 

.069159 

-.058445 

.026908 

-.005074 

.002776 

-.001627 

.031910 

-.011210 

.008927 

-.006836 

.492141 

-.070298 

.034939 

-.024007 

.647530 

-.172488 

.121519 

-.107639 

.036908 

-.003903 

.004824 

-.003375 

.043977 

-.019678 

.015509 

-.014228 

0.799339 

0.493783 

-0.089550 

0.049339 

1.282044 

1.145470 

-.310549 

.219544 

.169599 

.037952 

-.011949 

.006007 

.212792 

.083245 

-.038398 

.025243 

.699256 

.523828 

-.099723 

.058134 

1.138472 

1.219101 

-.345551 

.260447 

.157833 

.041704 

-.013166 

.008078 

.198610 

.091532 

-.042299 

.034055 

-0.313591 

0.651687 

0.575298 

-0.126091 

-.465823 

1.718204 

1.923065 

-.559573 

-.044602 

.150488 

.049347 

-.014691 

-.054326 

.340126 

.157843 

-.061693 

-.228783 

.633812 

.573549 

-.137821 

-.344915 

1.674943 

1.916360 

-.616234 

-.034985 

.148520 

.048803 

-.018600 

-.042759 

.335798 

.156076 

-.078385 

0.187897 

-0.159325 

0.576786 

0.562897. 

.276637 

-.400446 

2.109970 

2.432887 

.025479 

-.024868 

. .140460 

.042295 

.030950 

-.055302 

.458397 

.177006 

.117132 

-.136452 

.557359 

.664560 

.175138 

-.344108 

2.041363 

2.901646 

.017186 1 

-.021868 

.137234 

.063512 

.020954 

-.048664 

.447980 

.267043 , 

-0.069026 

0.050904 

-0.080137 

0.525349 

-.101352 

.127410 

-.283759 

2.458336 

-.009225 

.007667 

-.013645 

.134478 

-.011196 

.017031 

-.044065 

.573757 

-.031711 

.032307 

-.056459 

.432107 

-.047311 

.081140 

-.200039 

2.030260 

-.004593 

.005002 

-.009675 

.116059 

-.005596 

.011120 

-.031249 1 

.495663 
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TABLE VI - STATION NUMBERS 
n = 6 



\k 

JX 

1 

2 

3 

4 

5 

6 

M 


1. 172073 

0.391101 

-0.032371 

0.013323 

-0.010149 

0.008879 

N 


•638800 

.501856 

-.078978 

.046300 

-.045644. 

.048522 

P 


.210893 

.024685 

-.003894 

.001823 

-.001498 

.001143 

Q 


.051551 

.029221 

-.008595 

.005860 

-.006322 

.005949 

M* 


.676394 

.474177 

-.059129 

.026582 

-.018685 

.015649 

N» 


.441269 

.621067 

-.144759 

.092370 

-.083901 

.085903 

m 


.160476 

.034473 

-.007337 

.003631 

-.002691 

.002241 

[g 

1 

.041616 

.041021 

-.016206 

.011674 

-.011350 

.011694 

M 

2 

-1.066598 

0.853106 

0.468124 

-0.070505 

0.044513 

-0.035782 

N 


-.466718 

1.360105 

1.081893 

-.244062 

.199948 

-.195451 

P 


-.127634 

.176358 

.034411 

-.009203 

.006452 

-.004561 

Q 


-.026505 

.220969 

.075401 

-.029565 

.027216 

-.023740 

M* 


-.303948 

.731991 

.503742 

-.085145 

.049793 

-.038661 

N» 


-.164215 

1.186681 

1.169252 

-.294736 

.223322 

-.212149 

P* 


-.050692 

.162263 

.038896 

-.011103 

.007043 

-.005503 

m 


-.011384 

.203987 

.085309 

-.035665 

.029701 

-.028708 

M 

3 

1.150584 

-0.404200 

0.693794 

0.546418 

-0.133366 

0.089627 

N 


.489124 

-.597296 

1.822457 

1.822457 

-.597296 

.489124 

P 


.130870 

-.055956 

.156259 

.045293 

-.018484 

.011225 

Q 


.026719 

-.068060 

.352909 

.144808 

-.077925 

.058410 

M* 


.267118 

-.275585 

.662165 

.553477 

-.126314 

.083246 

N f 


.141177 

-.413819 

1.745286 

1.846431 

-.564890 

.456509 

P» 


.042839 

-.041308 

.152473 

.045980 

-.017100 

.011713 

Q* 


.009490 

-.050434 

.344556 

.146998 

-.072063 

.061097 

M 


-0.930965 

0.273028 

-0.194854 

0.592473 

0.624591 

-0.171416 

N 


-.390902 

.399897 

-.488124 

2.163786 

2.720209 

-.933437 

P 


-.103635 

.036020 

-.029594 

.142229 

.056698 

-.020561 

Q 

■ 

-.021011 

.043690 

-.065761 

.464052 

.238215 

-.106938 

M* 


-.210538 

.182630 

-,ie0822 

.598464 

.604424 

-.166643 

N* 


-.110263 

.271857 

-.454522 

2.185427 

2.628772 

-.912352 

Bl 


-.033225 

.026154 

-.028221 

.143453 

.053359 

-.022768 

m 


-.007320 

.031846 

-.062753 

.468015 

.224100 

-.118729 

M 

5 

0.581796 

-0.156399 

0.092907 

-0.111954 

0.537351 

0.599157 

N 


.242612 

-.228221 

.231501 

-.394888 

2.509279 

3.194001 

P 


.063996 

-.020218 

.013474 

-.018261 

.134578 

' .046991 

Q 


.012925 

-.024496 

.029896 

-.058920 

.574036 

.243745 

M* 


.120308 

-.097119 

.081568 

-.110987 

.535339 

.685021 

N* 


.062735 

-.144101 

.204039 

-.391731 

2.499695 

3.678582 

m 


.018841 

-.013674 

.012212 

-.016226 

.133988 

.066470 

hi 


.004140 

-.016634 

.027119 

-.058815 

.571521 

.345986 

M 

6 

-0.209220 

0.054246 

-0.030239 

0.031561 

-0.064035 

0.510543 

N 


-.086982 

.079042 

-.075223 

.110987 

-.291427 

2.902746 

P 


-.022893 

.006957 

-.004323 

.004969 

-.011243 

.132560 

Q 


-.004616 

.008425 

-.009587 

.016022 

-.047575 

.698254 

IV 


-.033141 

.025961 

-.020586 

.024431 

-.049677 

.425817 

N 1 


-.017248 

.038471 

-.051417 

.085977 

-.225999 

2.427820 

F» 


-.005173 

.003631 

-.003042 

.003877 

-.008676 

.115150 

V 


-.001135 

.004415 

-.006753 

.012502 

-.036708 

.606979 
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0.205449 

.083943 

.021819 

.004358 

.033023 

.016993 

.005053 

.001102 


0.137401 

.199169 

.017129 

.020712 

.086131 

.127049 

.011759 

.014281 


- 0.044613 

-.064008 

-.005533 

-.006689 

-.022307 

-.032878 

—.003033 

-.003682 


0.134585 

.334325 

.019010 

.042147 

.122950 

.306705 

.017983 

.039907 

- 0.064103 

-.158887 

-.008878 

-.019672 

-.058376 

-.145336 

-.008393 

-.010616 


0.020059 

.049675 

♦002756 

.006106 

.014642 

.036425 

.002091 

.004637 


0.629063 

2.291470 

.147488 

.480975 

.623577 

2.273022 

.147041 

.479557 


- 0.136596 
-.480675 
-.021700 
-.069981 
-.143020 
- . 503685 
-.022914 
-.073907 

0.054068 

.189539 

.008211 

.026453 

.057414 

.201491 

.008826 

.028441 


- 0.015652 

-.055505 

-.002373 

-.007644 

-.013579 

-.047602 

-.002060 

-.006638 


- 0.104856 

-.468955 

-.014225 

-.059956 

-.108454 

-.484684 

-.014528 

-.061218 


0.592219 

2.574726 

.051989 

.218352 

.584350 

2.538692 

.050504 

.212061 


0.551252 

2.570992 

.136202 

.580855 

.568583 

2.649577 

.139081 

.593022 

- 0.084474 

-.383331 

-.014243 

-.060230 

-.092005 

-.417397 

-.015450 

-.065330 


0.021382 

.096798 

.003481 

.014713 

.018939 

.085721 

.003074 

.012992 


- 0.027896 

-.153603 

-.004239 

-.022136 

-.031029 

-.170555 

-.004565 

-.023823 


0.081487 

.448232 

.012162 

.063490 

.073511 

.403671 

.010624 

.055439 


- 0.182666 

- 1.002357 

-.026098 

-.336173 

-.157315 

-.861833 

-.021773 

-.113560 


0.668960 

3.589325 

.063551 

.330701 

.634247 

3.397336 

.057851 

.300912 

0.507772 

2.883613 

.130010 

.684756 

.516450 

2.931074 

.131144 

.690659 


- 0.053078 

-.295079 

-.009551 

-.049981 

-.044219 

-.245636 

-.007854 

-.041094 


- 0.006513 

-.042160 

-.000863 

-.005357 

-.011099 

-.072080 

-.001614 

-.010037 


0.028701 

.185730 

.003780 

.023463 

.028988 

.188216 

.004199 

.026110 


- 0.077078 

-.498585 

-.010056 

-.062408 

-.063669 

-.413268 

-.009166 

-.056987 


0.144688 

.935220 

.018551 

.115110 

.115706 

.750618 

.016464 

.102355 


- 0.220971 

- 1.425675 

-.027154 

-.168426 

-.197770 

- 1.281188 

-.027322 

-.169818 


0.632193 

4.007039 

.051386 

.318021 

.704705 

4.491463 

.069349 

.430362 


0.498306 

3.334065 

.130932 

.820708 

.420133 

2.816717 

.114318 

.716958 



















































NACA TW 2185 


55 


TABLE VIII - STATION NUMBERS 
n « 8 



s 

1 

2 

3 

4 

5 

6 

7 

8 

u 

N 

P 

Q 

M» 

N* 

P' 

Q' 

1 

1.312192 

.695399 

.225483 

.054447 

.727233 

.467152 

.168111 

.043271 

0.364019 

.462793 

.021401 

.025255 

.444367 

.577362 

.030533 

.036245 

-0.021829 

-.052953 

-.002485 

-.005476 

-.043004 

-.104816 

-.005119 

-.011294 

0.006490 

.022453 

.000839 

.002694 

.015242 

.052799 

.002003 

.006433 

-0.00354 5 
-.015893 
-.000499 
-.002103 
-.008698 
-.038988 
-.001220 
-.005146 

0.003081 

.016958 

.000434 

.002422 

.007082 

.038933 

.001045 

.005456 

-0.003931 

-.025621 

-.000622 

-.003869 

-.007522 

-.048943 

-.001145 

-.007123 

0.005039 

.037691 

.000684 

.004929 

.008334 

.062492 

.001228 

.008863 

H 

N 

P 

Q 

M» 

N* 

P» 

Q' 

2 

-1.600390 
-.682210 
-.183160 
-.037524 
-.415706 ’ 

-.221089 
-.067464 
-.015018 

0.955663 

1.507998 

.188781 

.235968 

.797175 

1.282229 

.170868 

.214419 

0.428499 

.984086 

.029123 

.063696 

.468735 

1.082551 

.034087 

.074659 

-0.045077 

-.155335 

-.005550 

-.017807 

-.060778 

-.209723 

-.007611 

-.024428 

0.020351 
.091123 
.002808 
.011838 
.028715 
- .128559 
.003953 
.016664 

-0.016189 

-.089039 

-.002412 

-.012593 

-.021414 

-.117660 

-.003129 

-.016330 

0.019610 . 
.127790 
.003082 
.019183 
.021627 
-.140686 
.003272 
.020358 

-0.024337 

-.182002 

-.003289 

-.023705 

-.023284 

-.174560 

-.003421 

-.024694 

M 

N 

P 

Q 

M ' 
N* 
P* 
Q 1 

3 

2.337500 

.967933 

.254178 

.051181 

.462586 

.240602 

.072148 

.015838 

-0.630523 

-.923581 

-.083331 

-.101109 

-.388839 

-.579786 

-.056240 

-.068533 

0.780389 

2.036155 

.167794 

.378439 

.722423 

1.094495 

.160734 

.362854 

0.491625 

1.631308 

.037442 

.119546 

.512112 

1.702145 

.040069 

.127977 

-0.082648 

-.368956 

-.010886 

-.045861 

-.091426 

-.408100 

-.012016 

-.050617 

0.054514 

.299530 

.007975 

.041626 

.056607 

.310711 

.008120 

.042369 

-0.060500 

-.394067 

-.009418 

-.058605 

-.052469 

-.341167 

-.007864 

-.048918 

0.071477 

.534419 

.009602 

.069206 

.054058 

.405216 

.007911 

.057107 

u 

N 

P 

Q 

M' 

N» 

P* 

Q* 

4 

-2.630070 

-1.075644 

-.279850 

-.055949 

-.466670 

-.240468 

-.071591 

-.015627 

0.593584 

.061856 

.074709 

.090392 

.329868 

.487559 

.045534 

.055337 

-0.315718 

-.786292 

-.045646 

-.101283 

-.258293 

-.646305 

-.038018 

-.086223 

0.667195 
2.424356 
.152881 
.498324 
.650770 
2.367821 
. 150902 
.491981 

0.558819 

2.424356 

.046981 

.197214 

.561588 

2.436346 

.047159 

.197950 

-0.143453 

-.780292 

-.020073 

-.104712 

-.135711 

-.743106 

-.018614 

-.097075 

0.132430 

.861856 

.020253 

.126011 

.104892 

.681473 

,015449 

.096084 

-0.143916 

-1.075644 

-.019152 

-.138032 

-.100193 

-.750848 

-.014572 

-.105178 

u 

N 

P 

Q 

M* 

N» 

P' 

Q' 

5 

2.192995 

.890699 

.230544 

.045912 

.389428 

.199636 

.059201 

.012882 

-0.4 54018 
-.656778 
-.055975 
-.067648 
-.253524 
-.373346 
-.034316 
-.041657 

0.201514 
. .499217 
.027793 
.061577 
.165704 
.412443 
.023783 
.052749 

-0.175130 

-.614926 

-.027134 

-.087460 

-.170613 

-.599857 

-.026822 

-.086477 

0 . 584107 
2.718847 
.141097 
.601511 
.591124 
2,750891 
.142374 
.606915 

0.633389 

3.393592 

.058226 

.302906 

.614158 

3.287044 

.054969 

.285874 

-0.237005 

-1.539301 

-.034728 

-.215987 

-.191506 

-1.241724 

-.027027 

-.168026 

0.215982 

1.613222 

.028253 

.203600 

.155413 

1.164103 

.022340 

.161236 

V. 

N 

P 

Q 

M' 

N' 

P' 

Q. 

6 

-1.350286 
-.546005 
-.140865 
-.027983 
-.263362 
-.134575 
' -.039810 

-.008646 

0.265559 

.383370 

.032367 

.039093 

.163432 

.240192 

.021885 

.026552 

-0.108064 
-.267117 
-.014611 
-.032351 
-.098083 
- -.243600 
-.013810 
-.030611 

0.078152 

.273370 

.011560 

.037251 

.084195- 

.294907 

.012656 

.040765 

-0.102878 

-.466005 

-.016907 

-.071468 

-.115917 

-.525024 

-.019020 

-.080395 

0.520374 

2.952258 

.131525 

.692634 

.543778 

3.081687 

.135367 

.712714 

0.709714 

4.523995 

.070024 

.434540 

.662638 

4.216400 

.062213 

.385904 

-0.274441 

-2.046630 

-.034431 

-.248035 

-.230699 

-1.725908 

-.032195 

-.232314 

H 

N 

P 

Q 

N* 

P» 

Q‘ 

7 

0.654484 

.263837 

.067914 

.013468 

.130945 

.066783 

.019727 

.004280 

-0.124395 

-.179348 

-.015052 

-.018172 

-.079057 

-.116059 

-.010524 

-.012763 

0.048081 

.118705 

.006430 

.014232 

.045388 

.112606 

.006331 

.014029 

-0.031856 
; - .111252 
-.004626 
-.014089 
-.036054 
-.126105 
-.005328 
-.017155 

0.034801 

.157170 

.J05472 

.023114 

.041862 

.189085 

.006593 

.027846 

-0.066829 

-.370670 

-.011570 

-.060517 

-.077799 

-.431248 

-.013324 

-.069680 

0.484484 
3.239551 
.126311 
.791745 
.500211 
3.341811 
.128649 
. 806281 

0.662747 

4.867795 

.055545 

.399381 

.723527 

5.337321 

.072134 

.519815 

u 

N 

P 

Q 

M» 

H' 

P* 

Q* 

8 

-0.202414 

-.081470 

-.020947 

-.004150 

-.033107 

-.016869 

-.004979 

-.001080 

0.037821 
.054494 
.004560 
.005504 
.019722 
- .028937 

.002618 
.003174 

-0.014268 

-.035205 

-.001898 

-.004201 

-.011093 

-.027507 

-.001541 

-.003414 

0.009104 

.031774 

.001312 

.004221 

.008527 

.029809 

.001252 

.004029 

-0.009307 

-.041992 

-.001443 

-.006093 

-.009326 

-.042089 

-.001450 

-.006122 

0.015388 
.085181 
• .002573 
.013452 
.015058 
.083310 
.002493 
.013034 

-0.045167 

-.296574 

-.008294 

-.051711 

-.039785 

-.260993 

-.007173 

-.044710 

0.487926 

3.754721 

.129520 

.941457 

.414996 

3.198326 

.113558 

.825790 




TABLE IX - COMPARISON OF RESULTS 


56 


NACA TN 2185 







































Torsional deflection, 9/9, 
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NACA TN 2185 






Distance from root x 
Length of beam * l 

(d) Torsional vibrations; first mode, n = 3. 

Figure 2. - Comparison of theoretical mode shapes with 
mode shapes obtained by taking n stations along the 
beam for torsional vibrations. 



Torsional deflection, 0/0, 


NACA TN 2185 
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Distance from root x 
Length of beam * i 

(f) Torsional vibrations; third mode, n = 3. 

Figure 2L. - Concluded. Comparison of theoretical mode 
shapes with mode shapes obtained by taking n stations 
along the beam for torsional vibrations. 





Bending deflection, y/y 
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(c) Bending vibrations; second mode, n = 2. 



Distance from root x 
Length of beam * i 

(d) Bending vibrations; first mode, n = 3. 

Figure 3. - Comparison of theoretical mode shapes with 
mode shapes obtained by taking n stations along the 
beam for bending vibrations. 




Bending deflection, y/y. 
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Distance from root x 
Length of beam ' T 

(f) Bending vibrations; third mode, n = 3. 

Figure 3. - Concluded. Comparison of theoretical mode 
shapes with mode shapes obtained by taking n sta- 
tions along the beam for bending vibrations. 





Frequency ratio 
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Coupling coefficient, c 

Figure 4. - Variation of frequency ratio 0 with coupling 
coefficient c for several values of uncoupled frequency 
ratio Y. 
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